AD6227i>7 


AFWL-TR-65-132,  Vol.  I 


AFWL-TR 
65-132 
Vol.  I 


ABSORPTION  COEFFICIENTS  OF  HEATED 
AIR:  A  COMPILATION  TO  24,000°K 


Volume  I 

D.  R.  Churchill,  B.  H.  Armstrong, 
and  K.  G.  Mueller 

Lockheed  Palo  Alto  Research  Laboratory 
Lockheed  Missiles  &  Space  Company 

Palo  Alto,  California 
Contract  AF  29(601)-6320 

TECHNICAL  REPORT  NO.  AFWL-TR-65-132,  Vol.  I 


October  1965 


1/ 


1/ 


AIR  FORCE  WEAPONS  LABORATORY 
Research  and  Technology  Division 
Air  Force  Systems  Command 
Kirtland  Air  Force  Bose 
New  Mexico 


AFWL-TR-65-  1  32.  Vol.  I 


Research  and  Technoloj?y  Division 
AIR  FORCE  WEAPONS  LABORATORY 
Air  Force  Systems  Command 
Kirtland  Air  Force  Base 
New  Mexico 


When  U.  S.  Government  drawings,  specifications,  or  other  data  are  used  for 
emy  pui^jose  other  than  a  definitely  related  Government  procurement  operation, 
the  Government  thereby  incurs  no  responsibility  nor  any  obligation  whatsoever, 
and  the  fact  that  the  Government  may  have  formulated,  furnished,  or  in  any 
way  supplied  the  said  drawings,  specifications,  or  other  data,  is  not  to  be 
regarded  by  implication  or  otherwise,  as  in  any  manner  licensing  the  holder 
or  any  other  person  or  corporation,  or  conveying  any  rights  or  permission  to 
mtinufacture ,  use,  or  sell  any  patented  invention  that  may  in  any  way  be 
related  thereto. 

This  report  is  made  available  for  study  with  the  understanding  that 
proprietary  interests  in  and  relating  thereto  will  not  be  impaired.  In 
case  of  apparent  conflict  or  any  other  questions  between  the  Government's 
rights  and  those  of  others,  notify  the  Judge  Advocate,  Air  Force  Systems 
Command,  Andrews  Air  Force  Bass,  Washington,  D.  C.  20331. 

DDC  release  to  OTS  is  authorized. 


AFWL-TR-65-132.  Vol.  I 


ABSORPTION  COEFFICIENTS  OF  HEATED 
AIR:  A  COMPILATION  TO  24,000OK 


Final  Report 
Volume  1 

O.  R.  Churchill*  B.  H.  Armstrong, 
and  K.  G.  Mueller 


Lockheed  Palo  Alto  Research  Laboratory 
Lockheed  Missiles  It  Space  Company 
Palo  Alto*  California 
Contract  AF  29(601). 6320 


TECHNICAL  REPORT  NO.  AFWL-TR-65-  132*  Vol.  I 


) 


AFWL-TR-65-132.  Vol.  I 


FOREWORD 


This  report  was  prepared  by  members  of  the  Physical  Sciences  Laboratory 
of  Lockheed  Missiles  and  Space  Company,  Palo  Alto,  California,  under  Contract 
AF  29(601)-6  320 .  The  work  was  performed  under  Program  E  lement  7. 60,  06. 01. D, 
Project  5710,  Sub  task  07.003,  and  was  funded  by  the  Defense  Atomic  Support 
Agency  (DASA).  This  effort  was  carried  out  at  the  Palo  Alto  Research  Labora¬ 
tories  during  the  period  from  1  April  1964  to  30  July  196  5.  The  report  was 
submitted  on  20  September  1965  by  the  AFWL  Project  Officer,  Lt  Fred  J.  Reule 
( WLRTH). 
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ABSTRACT 

An  extensive  tabulation  of  absorption  coefficients  of  heated  air  has  been 
carried  out  with  the  aid  of  several  large  digital  computer  programs.  These 
tables  supersede  the  previous  ones  computed  at  Lockheed  by  Meyerott, 
Sokoloff,  and  Nicholls  in  1959.  Tables  are  presented  for  temperatures  at 
every  1000°K  between  1000°K  and  24,  000°K  and  for  eight  densities  between 
ten  times  normal  atmospheric  to  10“^  times  normal.  Absorption  coefficients 
are  individually  listed  for  each  of  14  contributing  absorbers.  The  photon 
energy  range  is  0.6>10.  7  eV  with  a  partition  interval  of  0.  1  eV. 
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INTRODUCTION 


The  taUes  of  absorption  coefficients  of  heated  air  presented  in  this  report  are  intended 
to  supersede  the  previous  compilation  of  Meyerott»  Sokoloff,  and  Nicholls  (Ref.  1). 

Their  format  for  presentation  of  the  absorption  coefficienU  will  be  followed  here,  since 
it  allows  the  user  to  convenienUy  select  the  important  contributors  from  the  listed 
individual  absorbers  at  any  temperature,  density,  and  photon  energy.  Differences 
between  the  two  calculations  are  extensive,  and  reflect  the  advances  made  in  the  field 
of  opacity  since  the  MSN  calculations  were  performed  in  1959.  Improved  com¬ 

puter  techniques  have  been  used  to  theoretically  reconstruct  the  line  spectra  for  the 
molecular  band  systems  which  are  included  in  the  present  tables,  where  the  previous 
calculation  had  lumped  the  entire  band  contribution  at  band-head  energy.  Recent 
experimental  results  and  analyses  have  made  available  improved  f  numbers  and  Franck- 
Condon  factors.  Atomic  lAotoionization  calculations  have  been  advanced  significantly 
since  1959  by  the  work  of  a  number  of  investigators,  allowing  us  to  greatly  improve 
over  the  crude  version  of  the  hydrogenic  model  of  photoionization  employed  by  MSN. 

Since  the  present  tables  are  the  result  of  a  series  of  large  digital  computer  programs, 
it  has  been  relatively  simple  to  extend  the  tables  to  24, 000*  K,  thus  providing  much 
needed  information  in  regions  where  little  has  been  available,  and  to  refine  the  energy 
partition. 

The  contributors  to  the  heated  air  absorption  which  have  been  included  in  our  calcula¬ 
tion  are  listed  in  Table  I . 
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Table  I.  Contributors  to  the  Absorption  Coefficient  of  Heated  Air 


1.  Schumann-Runge  bands  of 

2.  First  positive  bands  of  N2 

3.  Second  positive  bands  of  N2 

4.  First  negative  bands  of 

5.  Beta  bands  of  NO 

6.  Gamma  bands  of  NO 

7.  Vibration- rotation  bands  of  NO 

8.  Birge-Hopfield  No.  1  bands  of  N^ 


9.  Schumann-Runge  continuum  of  O2 

10.  NO2 

11.  0~  photodetachment 

12.  Free-free  in  the  presence  of  ions 


13.  Atomic  nitrogen  photoionization 

14.  Atomic  oxygen  photoionization 


These  may  be  separated  into  three  groups  according  to  their  method  of  treatment  in 
the  calculaticm.  Contributors  one  through  ei^t  arise  from  bound-bound  transitions  in 
diatomic  molecules  and  were  treated  throu^  the  extensive  SACHA  (Refs.  2  —  4)  digital 
computer  programs.  Contributors  13  and  14  were  con\puted  with  the  large  PIC  compu¬ 
ter  progpram  (Ref.  5).  The  rest  were  entered  into  the  calculation  in  a  semi -empirical 
fashion  and  might  be  called  continuum  and  pseudo- continuum  processes,  even  though 
discrete  transitions  may  have  been  included  for  NO2,  (Ref.  6).  A  special  computer 
code  was  written  to  collate  the  results  into  the  desired  format  and  print  out  the  tables. 
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Theoretical  construction  of  the  contributions  from  bands  and  continua  of  diatomic 
molecules  is  discussed  in  Section  n,  absorption  arising  from  photoionization  of  atomic 
nitrogen  and  oxygen  is  discussed  in  Section  m,  and  the  remaining  processes  are 
covered  in  Section  IV .  Comments  on  the  results  are  contained  in  Section  V.  The  bulk 
of  the  report  is  contained  in  the  form  of  tables  of  absorption  coefficients;  these  appear 
in  Volume  II .  The  equilibrium  code  which  was  used  in  these  calculations  is  discussed 
in  Appendix  A. 
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Section  II 

ABSORPTION  BY  DIATOMIC  MOLECULES 
2. 1  LINE  ABSORPTION 

llie  line  absorption  coefficient  for  a  radiative  transition  between  two  nondegenerate 
levels  L  (lower)  and  U  (upper)  may  be  written  (Ref.  7) 

'  N^»n,hi>n,F(i.)  (2.1) 

where  ie  the  Bohr  frequency  of  the  line  In  cm"*,  la  the  particle  number  denalty 
(number  of  particles  per  unit  volume)  in  the  lower  level,  Bj^^  is  the  Einstein  coefficient 
for  absorption,  and  F(v)  is  a  line  shape  factor  normalized  such  that  /f(v)  dv  =  1. 

To  determine  the  quantities  on  the  rlght'^and  side  of  Eq.  (2, 1),  the  energy  levels, 
selectian  rules,  and  wave  functions  (or  transition  probabilities)  must  be  found  by 
application  of  the  laws  of  quantum  mechanics  to  the  system  under  investigation,  which 
in  this  case  is  a  diatomic  molecule.  The  theoretical  approach  is  to  effect  a  separation 
of  Schrbdinger's  equaUon  by  assuming  there  to  be  no  interaction  between  the  modes  of 
nuclear  vibraUon.  rotaUon,  and  electronic  moUon.  This  method  is  known  as  the 
Bom-OppenheimerapproximaUon.  Once  the  necessary  theoreUcal  expressions  are 
obtained,  they  may  be  corrected  for  small  interactions. 

Theoretical  investigations  on  diatomic  molecules  are  less  complicated  then  those  on 
polyatomic  nK>leculcs  because  of  the  added  vibrational  modes,  rotational  asymmetries, 
and  the  larger  number  of  possible  electronic  states  generally  associatcnl  with  the  latter. 
However,  purely  theoretical  approaches  yield  much  less  satisfactory  results  for  dia¬ 
tomic  molecules  than  for  atoms,  particularly  w'hen  accurate  wave  functions  are  desired. 
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The  electronic,  vibrational,  and  rotaticmal  energies  of  a  diatomic  molecule  are 
quantized  (we  are  not  concerned  with  the  unquantized  translational  energy  in  the  radia¬ 
tion  problem).  Of  the  set  of  quantum  numbers  which  specify  a  given  level,  the  msg- 
netic  quantum  nunriber  M  (the  Z  compcment  of  the  total  angula**  momentum  quantum 
number  J)  may  be  omitted  from  the  expression  for  the  absorption  coefficient  which  is 
applicable  to  a  gas  of  randomly  oriented  absorbers  since  the  magnetic  quantum  number 
specifies  the  orientation  in  qpaoe  and  is  therefore  eliminated  by  the  usual  procedure 
of  averaging  over  initial  sub-levels  and  summing  over  the  final  ones.  In  addition, 
the  spin  quantum  number  S  enters  in  a  nonessential  way  since  spin  is  unaffected  in  the 
electric  d4>ole  transitions  we  are  c<msidering  here. 

The  Einstein  Coefficients  Electronic  transitions  in  diatomic  molecules  will,  in 

accord  with  spectroscopic  convention  (Ref.  8),  be  labeled  by  n*,v',  J*  -  n**,v*',J*', 
where  n  denotes  the  electronic  state  and  v  and  J  are  the  vibrational  and  rotational 
cpiantum  numbers,  respectively.  Double  primes  refer  to  lower  levels  and  single 
primes  to  upper  levels.  As  applied  to  transitions  in  diatomic  molecules,  the  Einstein 
coefficient  for  absorption  is  given  by  (Ref.  8) 
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2  IrM'.M" 


n”,v**,J”,n»,v»,J» 


8ir  M».M” 

3h2^  2J**  ^  1 


(2.2) 


udiere 


^  J^n*’,v”,J**,M**  “  *n»,v*,J»,M»  ® 


(2.3) 


is  the  matrix  element  of  the  electric  d4x>le  moment  R  »  R  R  of  the  electrons 

e  nuc. 

and  imclei.  M*  and  M"  label  die  spatially  degenerate  components  of  the  upper  and 
lower  rotational  levels  and  the  summation  is  over  all  possible  oonibinatians  of  M* 
and  M"  . 
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We  assume  validity  of  the  Born-Oi^nheimer  f^sproximation  in  the  form  wherein  the 
molecular  wave  function  may  be  e^q^ressed  as  a  product  of  electronic,  vibrational, 
and  rotational  wave  functions,  viz: 


♦  ♦  ♦ 
e  V 


rot 


(2.4) 


with  this  approximation, the  expression  for  the  matrix  element  becomes 


.M 


e  V  rot' 


=  /(*;'*rot>  ««<■•>  <n*rot>  ^’’vrot 


(2.5) 


where  R^(r)  =  /♦”  R  dr^  is  the  electronic  electric-dipole  transition  moment, 
and  is  assumed  to  vary  slowly  with  the  nuclear  separation  coordinate,  r  .  Since 

the  electronic  wave  functions  are  orthogonal,  the  nuclear  part  of  R  does  not  contribute 
to  R^(r). 


The  sum  in  Eq.  (2.2)  may  now  be  expressed  in  the  form  (cf,  Penner  (Ref.  9,  p.  129)J 


2  |rM’.M'12 

M'.M" 


(v"|Re(r)!v)2Sj,._j, 


(2.6) 


where  j,  is  the  so-called  Hbnl-London  factor,  which  contains  all  the  angular 
dependence.  Following  Nicholls  and  Stewart  (Ref.  10),  we  adopt  the  term  band  strength 
for  the  other  factor,  as  suggested  by  Chamberlain  (Ref.  11),  with  the  definition 

(2.7) 
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According  to  Fraser  (Ref.  12),  it  is  possible  to  rewrite  this  as 

under  many  circumstances,  and  for  many  band  systems.  In  Eq.  (2. 8), 

q(v*.v’’)  £ 

is  the  Franck- Ccmdon  factor  for  the  band  and 


(2.  8) 


(2.9) 


,r*  ,,  dr 

-  ^  J  V*  v” 

v’.v"  -  dr 


(2. 10) 


is  the  corresponding  r- centroid  (Ref.  13). 

Accurate  calculations  of  the  electronic  transition  moments  are  not  available  for  the 
molecules  of  interest  to  this  study.  In  general,  lacking  these  calculations,  recourse 
must  be  made  to  experiment  for  the  determination  of  the  "absolute"  values  of  these 
quantities.  Experimental  band  strength  information  is  not  available  for  most  of  the 
bands  present  in  heated  air,  however,  and  the  desired  S  ,  ,,'s  must  be  estimated. 

These  estimates  are  usually  obtained  in  the  following  manner. 

From  the  available  measurements  for  one  or  more  bands  (v'.v")  of  a  band  system, 

and  with  the  aid  of  calculated  Franck-Condon  factors,  the  R  (r  „)*s  are  calculated 

e  V  * ,  V 

by  means  of  Eq.  (2.8).  Using  these  values  of  R  ,  a  smooth  curve  is  drawn  from 
which  the  desired  values  of  are  obtained  for  the  various  bands.  Multiplica¬ 

tion  by  the  appropriate  Franck-Condon  factors  then  yields  a  corresponding  set  of  band 
strengths,  .  This  procedure  is  practical  because  R^(r)  is  usually  a  relatively 

slowly  varying  function  of  r  .  (Franck-Condon  factors,  on  the  other  hand,  vary 
over  several  orders  of  magnitude  for  each  band  system  considered  here. )  In  a  recent 
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review  article,  Nicholls  (Ref.  14)  gave  an  extensive  discussion  of  the  basic  theory 
outlined  above  and  presented  tables  of  band  streng^s  for  many  band  systems  of  air 
molecules.  We  have  drawn  heavily  upon  this  article  for  the  present  study. 

The  procedure  for  obtaining  the  matrix  elements  fEq.  (2.3))  for  vibration- rotation 
transitions  in  diatomic  molecules  differs  from  that  discussed  above,  and  is  fully 
discussed  for  the  one  vibration-rotation  band  system  included  in  the  present  calculation 
by  Churchill  and  Hagstrom  (Ref.  4). 

The  number  densities  or  occupaticm  numbers  Nj^.  The  particle  number  density  in  the 
lower  level  (often  called  the  occupation  number)  is  given  by  the  usual  expression 
(Ref.  15) 


N 


n*',v",J” 


N(2J*'  l)a>j  ejqp 

Q 


[  ^"kT  '^  ] 


(2. 11) 


where  N  is  the  total  number  density  for  all  levels,  E^„  is  the  energy  of  the 

n  f  V  9  u 

level  and  o;^  its  nuclear  spin  statistical  weight.  Q  is  the  total  partition  function  and 
may  be  written  as  the  sum  of  contributions  from  all  the  electronic  states: 


(2. 12) 


Taking  into  account  electronic  (orbital  spin),  vibrational,  rotational,  and  nuclear 
(spin  only)  degrees  of  freedom,  we  may  write 


where 


^electronic  ^vib-rot  Nuclear 


^electronic  [  kT  | 


(2. 13) 


(2. 14) 
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I 


^vlb- 


rot 


=  2  exp  [ 
v=0 


-G 


kT  J  '^rot 


(2. 15) 


Q 


(V) 

rot 


•e 

■/ 


(2J  +  1)  exp 


^-Fy(J)hc| 


dJ 


kT 


^  <2-16) 


(21  +  1)  (2L  +  1) 

O  ,  =  - S - US - 

niuclear  o 


(2. 17) 


In  these  formulas 


is  the  energy  of  tho  v 
oo 


ft  — 


=  0  to  the  v'  =  0  transition  of  the 


electronic  states  in  question.  G  (v)  is  the  vibrational  term  referred  to  the  lowest 


vibrational  level,  F^(J)  is  the  rotational  term  for  the  v 


th 


vibrational  level  with  B 


its  corresponding  rotaUonal  constant,  and  S  is  the  total  electron  spin,  is  the 
statisUcal  weight  for  orbital  angular  momentum  A  about  the  intemuclear  axis  and 
takes  the  value  1  if  A  =  0  and  value  2  if  A  0.  and  Ij^  are  the  nuclear  spins 
of  the  two  nuclei  and  a  is  a  symmetry  number  having  the  value  2  for  homonuclear 
molecules  and  1  for  heteronuclear  systems. 


We  have  made  use  of  an  approximate  formula  for  Qyj^-rot  ^ore  convenient 

for  calculation  and  yet  accuraie  enough  for  our  problem.  This  formula  was  developed 
by  Bethe  (Ref.  16)  and  later  corrected  by  Brinkley,  et  si.  (Ref.  17). 


^vib- 


rot 


1  -  exp 


( 


-1.4388a; 


4388B. 


(1  +  VT) 


(2. 18) 


where  y,  Bethe 's  correction  factor  for  anharmonicity  and  non-rigidity,  has  the  value 


y  = 


1.4388a; 


_2_o  +  o  ^  __o\ 

a;  B^  a;  / 

o  O  Of 


(2.  19) 
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Here  and  are  the  vibrational  and  rotational  conatanta  «.€  the  v  =  0  vibrational 
level,  ia  the  firat  anharmonic  conatant  and  ia  the  interaction  conatant  for 

coupling  between  rotation  and  vibration. 


Returning  now  to  Eq.  (2, 11),  we  multiply  nunterator  and  denominator  by  Q  aa 
given  by  Eqa.  (2. 13),  (2. 14),  and  (2. 17)  and  obtain  after  a  little  algebra 


^n*V”J"  ■  ^ 


total  Q, 


%>>  *  1)  e«P  I  -  {Go<v")  *  Fy.4J">) 


total 


<2S  ^^‘‘^A^uclear^ib  -rot 


(2.  20) 


where  uae  has  been  made  of  the  relation 


®n*V*’J”  "  ^‘'oo  rvt.(J”))  he  .  (2.21) 

Now  let  =  Si”^^total  **®"®*®  *•*«  fractional  population  of  the  n"***  elecrtronic 
state  and  introduce  the  notation 


8ir^L^  Wj 

"n’n"  VV-.J’-*"  °  3hc  (2S  +  n'n" .  J’J"  J-  22, 

and 

<2-23> 

where  =  2.6875  x  lo^®  particlea/cm^  ia  Loachmidt 'a  number. 

We  now  obtain  the  desired  form  of  Etj.  (2. 1)  by  substitution,  using  Eqa.  (2.20),  (2.22), 
and  (2.  23); 


M(i') 


L^o^ib-rotJ 


v'v'*. 


rF(i#) 


(2.  24) 
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Here  the  quantities  H  and  E  are  characteristic  of  the  isolated  molecule;  i.e. , 

do  not  depend  on  the  temperature  and  density.  Since  the  bracketed  term,  in  (2.24)  is 

-2 

dimensionless,  H  has  dimensions  cm  ,  and  F(p')  has  dimensions  cm,  the  absoi^ption 
coefficient  p(r)  has  the  r-'quired  dimension  cm  ^ . 

Tlie  foregoing  equations  have  been  incorporated  into  the  digital  computer  programs 
which  reconstruct  the  absorption  spectra  of  the  important  band  systems  of  heated  air. 
These  band  systems  will  be  discussed  a  little  further  on. 

Species  concentrations  and  Boltzmann  factors  for  the  electronic  states  of  diatomic 
air  molecules  are  computed  by  a  separate  program  known  as  the  Equilibrium  code. 

The  results  of  this  computation,  which  is  discussed  in  Appendix  B,  are  provided  to  the 
other  programs  on  a  magnetic  tape  for  digital  computer  use. 

The  average  absorption  coefficient  and  the  line  shape  factor  F(r).  The  absorption 
coefficient, as  required  in  the  strict  development  of  the  equation  of  radiative  transfer, 
is  a  continuous  function  of  frequency.  However, it  is  rarely  possible  to  make  use  of 
such  monochromatic  values  of  the  absorption  coefficient  in  practical  problems  concerning 
the  flow  of  radiative  energy  from  one  location  to  another  because  of  the  mathematical 
complexity  of  such  problems.  Their  practical  solution  is  usually  undertaken  \^ath  the 
aid  of  appropriately  chosen  "average"  absorption  coefficients  which  are  assumed  con¬ 
stant  over  some  suitable  frequency  interval.  In  addition,  the  shape  of  molecular  lines, 
as  represented  by  the  line  shape  factor  F(r),  are  not  yet  well  Imown  as  functions  of 
temperature  and  density.  A  number  of  specific  line  broadening  mechanisms  have  been 
studied  (Ref.  18)  both  experimentally  and  theoretically.  However,  the  uncertainties 
of  the  general  problem  are  such  as  to  make  it  inadvisable  to  include  such  partial  infor¬ 
mation  in  an  absorption  coefficient  study  which  covers  such  a  broad  range  of  tempera¬ 
tures  and  densities  as  does  the  present  one.  For  the  reasons  outlined  above  we  make 
use  of  a  frequency  averaging  procedure  in  our  calculations  so  that  we  do  not  present 
monochromatic  absorption  coefficients.  Even  with  this  enormous  simplification , the 
amount  of  data  to  be  computed  and  presented  is  very  large  and  requires  extensive  use 
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of  large  digital  computer  facilities.  Again,  due  to  the  complexity  of  the  radiative 
transfer  equation,  it  is  not  possible  to  specify  a  priori  the  nature  and  magnitude  of 
the  errors  which  will  arise  from  our  frequency  averaging  procedures.  Such  errors 
will,  of  course,  depend  on  the  particular  application,  and  the  reader  is  cautioned  to 
bear  in  mind  the  limitations  inherent  in  the  use  of  our  results,  particularly  in  connec¬ 
tion  with  problems  involving  large  optical  depths. 

The  averaging  procedure.  Suppose  o?  denotes  a  set  of  numbers  which  completely 
specifies  the  transition  from  which  a  given  line  arises.  We  are  concerned  with  a 
large  but  finite  collection  of  such  sets  of  numbers  and  so  we  can  order  these  arbitrarily. 
Labeling  the  line  L-U  with  the  vector  a  ,  Eq.  (2.11)  may  be  written 

M  (t')  =  N  B  hP  F  (i^)  (2.25) 

O’'  '  a  a  a  a'  '  \  / 


The  reason  for  the  change  in  notation  becomes  clear  when  we  define  an  average  absorp¬ 
tion  coefficient  which  pertains  to  a  frequency  interval  Av : 


(2.  26) 


and, since 


o 


where  p  is  the  contribution  from  continuum  transitions,  then 
c 

Au  " 


Au 


I 


N  B  h 

a  O' 


■'  f 

a  J 


F  (u)  dv  + 
o'  '  c 


Au 


(2.27) 


(2.28) 
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We  now  assume  that  lines  whose  Bohr  frequencies  lie  outside  the  frequency  interval 
Av  do  not  contribute  to  the  average  absorption  coefficient  for  that  interval.  The  error 
introduced  by  this  assumption*  is  reduced  by  taking  Av  large  enough  so  that  many 
lines  are  contained  therein.  Furthermore,  the  average  defined  in  Eq.  (2. 26)  is 
generally  used  only  in  calculations  which  apply  to  optically  thin  materials  where  the 
lines  remain  relatively  narrow,  and  so  do  not  contribute  appreciably  to  frequencies 
whi<^  differ  from  the  B<rfir  frequency. 

It  is  then  reasonable  to  assume  that 


/  ^o^*'^***'  *  /  =  1  (2.29) 

Ap  o 

and  Eq,  (2. 28)  takes  on  the  simple  form 

^  Z  ♦  5^  (2,30) 

where  the  sum  is  taken  over  those  lines  whose  maxims  fall  within  Av  . 

The  sum  appearing  in  the  preceding  expression  may  be  partitiooed  into  sub-sums 
acoording  to  molecular  species,  particular  band  system,  etc.  In  this  report  we  are 
concerned  with  oontributioos  from  individual  band  systems  and  the  appropriate  parti¬ 
tion  has  been  made. 


*ln  the  sense  of  a  direct  average,  not  a  harmonic  average  such  as,  for  example,  the 
Rosseland  mean.  For  this  latter  type  of  average  our  procedure  actually  increases 
the  error. 
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2.2  THE  CONTRIBUTING  BAND  SYSTEMS 
2.2. 1  Band  ^atema  Included  in  the  Calculation 

The  molecular  hand  ayatema  which  have  been  included  in  the  calculation  are  liated  in 
Table  II.  Detalla  of  the  theoretical  reconatruction  of  all  of  theae  band  ayatema  except 
one  have  been  given  in  previoua  publicationa  (Refa.  2-4)  and  therefore  wUl  not  be 
repeated  here.  The  exception  ia  the  Birge-Hopfield  No.  1  ayatem  of  N2  which  will 
be  diaouaaed  next.  It  ahould  be  noted*  however*  that  one  of  the  above-cited  refeiwnoea 
(Ref.  3)  givea  an  eiqireaaion  for  the  line  intenaity  factor*  to  be  ^n>lied  to  the  N^  firat 
poaitive  ayatem.  which  containa  HBal- London  tectora  that  do  not  apply  to  n»et  ayatem 
(i. e.*  they  are  incorrect).  The  correct  Hbnl-London  factora  were  actually  uaed  in 
the  calculationa  and  are  given  ia  Ref.  2. 


Table  n.  Molecular  Band  ^  atema  Included  in  the  CaloulaUon 
of  the  Abaorption  Coefficient  of  Heated  Air 


Band 

Syatem 

Number 

Molecular  Speciea 
and  Syatem  Name 

Tranaition  in 
Spectroscopic 
Natation 

Total  Number  of  Linea 
Included  in  the  Calculation 
for  Ihia  Band  ^atem 

1 

O2  Schumann- Runge 

.  X’S* 

4*611 

2 

N2  Oral  poaitive 

B®n^  .  A  *2* 

48*785 

3 

N2  second  poaitive 

c*n„  -  B*n^ 

20*369 

4 

N^  firat  negative 

B*r;  - 

3*216 

5 

NO  beta 

B  -  X 

18*518 

6 

NO  gamma 

A  -  X 

31*429 

7 

NO  vibration- rotation 

X  -  X  *11 

25*610 

8 

^2  Birge-Hopfield  No.  1 

b*‘r*  -  x*s* 

38,983 

14 


f 


The  spectral  atlases  for  all  of  these  systems  have  been  merged  into  a  master  air 
atlas  to  be  used  in  this  and  other  calculations.  This  atlas  contains  the  following 
information  for  each  line  to  be  included  in  the  air  computation:  the  frequency,  H 
factor  (described  previously),  and  identification  vector.  The  identification  vector 
permiU  lines  to  be  selected  according  to  their  band  system  so  that  the  proper  tempera¬ 
ture-density  dependent  factor  may  be  applied  to  account  for  the  parUcle  populaUon 
density  in  Ute  lower  state  of  the  transition. 

The  band  strengths  urhich  were  used  in  the  reconstruction  of  the  first  6  ^sterns  listed 
in  Table  U  were  taken  directly  from  a  compilation  by  NichoUs  (Ref.  14).  Otherwise, 
the  qpectrosoopic  constanU  used  in  the  reconst:  ctlon  of  the  first  seven  listed  systems 
renudn  the  same  as  for  previous  calculations  (Refs.  2-4). 

2.2.2  The  Birge-Hopfleld  No.  1  Band  System  of  N^ 

The  vacuum  ultraviolet  spectrum  of  N^  to  very  complicated,  with  a  large  number  of 
band  systems  having  been  identified  (Ref.  it).  flk>eotrosoopic  constaato  are  available 
for  some  of  these,  indudiag  the  Birge-HopOeld  No.  1  beads  (Ref.  20)  which  arise 
from  transitkms  between  the  ground  electronic  stale  X  and  the  eacited  state 
b*  .  Sinoe  these  bands  are  relatively  strong  to  absorption  and  extend  to  longer 
wavelengths  (*«  1450  A)  than  most  of  the  others,  they  were  includsd  to  the  present 
calculation.  The  Birge-Hopfleld  No.  2  band  system,  which  lies  to  the  same  frsqtteac>' 
region,  was  not  includsd.  Farther  reaMrks  on  this  ondaslon  will  be  reserved  for 
Section  V. 

Since  both  states  invohrsd  to  the  Birge-HopOeld  No.  1  transitions  are  *s,  the 
structure  of  the  band  syalem  to  simple.  Each  band  conaisU  of  a  stogie  R  md  a 
single  P  branch  (since  AJ  -  O  to  forbidden  for  Z  -  Z  transitions,  there  is  no  Q 
branch). 


15 


Line  frequencies  for  these  bends  may  be  ceiculeted  by  use  of  the  relation 


Vrt'WV.J’J’*  '  •'oo  *  <2.31) 

*00  ^•**»®*®™*  for  the  v”  -  0  v*  =  O  rotationless  transition,  G  C*)  is 
the  vibrational  term  of  the  level  where  the  first  vibrational  level  is  v  «  O  .  mk) 
is  the  rotational  term  for  the  v^  level.  For  the  R  branch  terms,  J*  >  J  1: 
for  the  P  branch  terms,  J*  J**  •  1  . 

Rotational  terms  are  computed  for  both  states  from  the  expression  (Ref.  8) 

‘  ■*'*<•'  ♦  *>  -  ♦  »*  (2.32) 

Observed  and  calculated  values  of  the  rotational  oonstanU  and  are  (iven  by 
Lofthus  (Ref.  20)  for  the  X  sUte  up  to  v**  -  18 .  He  also  gives  G  (v)  iq>  to 
v”  -  20  .  ® 


Ws  have  used  the  formula  of  Wilkinsen  as  givsn  by  Lofthus  (Ref.  20)  to  obtain  values 
of  for  V**  up  to  27  for  the  X  state,  and  we  have  used  the  formula  of  Lofthus 
(Ref.  20)  to  obtain  values  of  G^(v) .  also  to  v**  -  27.  was  assumed  oonaUmt  (at 
the  observed  value  of  •  x  io-«  cm-1  for  v”  -  18)  for  all  v”  e  18  . 


Table  10  lisU  the  apeetroacoplc  constants  for  the  X  electmnic  state  which  were 
used  in  our  caloulaUon.  was  taken  to  be  103, 072. 1  cm*^  as  givsn  by  Lofthus 

(Ref.  20).  The  maximum  aUowabIs  value  of  J  was  set  to  100  for  our  ^oulatloa. 


*The  word  **term**  is  used  in  molecular 


to 


in  the  dimension  cm*l 


I 

i 


Table  in.  Spectroscopic  Constants  for  the  X  State  of 


V 

Go(vl 

(cm”l) 

By 

(cra“l) 

Dv 

(cm'l) 

0 

00.00 

1.9898 

6.4  X  lo"** 

I 

2329.66 

1.9720 

5.0  X  10"® 

2 

4630.83 

1.9548 

5.0  X  lO’®* 

3 

6903.44 

1.9364 

5.5  X  10’®* 

4 

9147.41 

1.9186 

5.5  X  10’®* 

5 

11362.67 

1.9006 

5.5  X  10’®* 

6 

13549. 14 

1.8829 

5.0  X  10"®* 

7 

15706.76 

1.8651 

5.0  X  lO"®* 

8 

17835.44 

1.8473 

5.0  X  10"®* 

9 

19935. 11 

1.8295 

5.0  X  lO"®* 

10 

22005.70 

1.8117 

5.0  X  10"®* 

11 

24047.13 

1.7939 

5.3  X  10“®* 

12 

26059.33 

1.7761 

5.5  X  10"®* 

13 

280U.23 

1.7583 

5.5  X  10"®* 

14 

29995.74 

1.7405 

6.0  X  lo"®* 

15 

31919.80 

1.7227 

6.0  X  10’®* 

16 

33814.33 

1.7048 

6.0  X  10’® 

17 

35679.26 

1.6870 

6.0  X  lO"® 

18 

37514. 51 

1.6892 

6.0  X  10“® 

19 

39320.01 

1.6514 

6.0  X  10’® 

20 

41096.68 

1.6338 

6.0  X  10’® 

21 

42841.45 

1.6158 

6.0  X  10”® 

22 

44557.25 

1.5980 

6.0  X  10’® 

23 

46242.99 

1.5802 

6.0  X  10"® 

24 

47898.61 

1.5624 

6.0  X  10’® 

.5 

49524.03 

1.5446 

6.0  X  10”® 

26 

51119. 17 

1.5267 

6.0  X  10’® 

27 

52683.97 

1.5089 

6.0  X  10’® 

^Observed  value. 
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For  the  b*  ^2^  upper  electronic  state,  the  spectroscopic  constants  were  calculated 
from  formulas  due  to  Lofthus  (Ref.  20).  which  give  good  agreement  with  the  available 
experimental  data.  Of  the  eight  observed  vibrational  levels  in  this  state,  seven  have 
been  included  in  our  calculation. 

The  aingle  omission  is  the  highest  (v**  =  7)  observed  level;  no  Franck-O>ndon  factors 
were  available  for  this  level.  Table  IV  lists  the  constants  emplo3red  in  the  calculations 
for  this  state. 


Table  IV.  Spectroscopic  Constants  for  the  b*  State  of  N2 


V 

Oo(v) 

(cm“f ) 

Bv 

(cra'l) 

Dv 

(cm~l> 

0 

0.0 

1.1515 

1.104  X  10’* 

1 

742.0 

1.146 

1.137  10"* 

2 

1474.4 

1.142 

1.169  X  10"* 

3 

2196.7 

1.137 

1.202  K  10"* 

4 

2909.6 

1.132 

1.234  X  10"* 

5 

3613.2 

1.128 

1.267  X  10"* 

6 

4307.0 

1.123 

1.299  10"* 

We  now  turn  to  a  discussion  of  the  line  iatensities  associated  with  these  bands. 

Since  nitrogen  nuclei  obey  Bose  staUeUcs  (1  -  1).  the  symmetric  levels  wlU  have  the 
higher  sUtistlcal  weight  (Ref.  S) 

Wj(s)  -  (21  1)  (1  4^  1)  -  6 

and  the  anUs>'mmetrlc  levels  will  have  the  lower  statistical  weight 

Wj(a)  •'  (21  ^  1)  I  -  3 
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For  the  state,  rotational  levels  belonging  to  even  J  are  symmetric  and  those 
belonging  to  odd  J  are  antis3nnmetric.  Consequently,  there  is  an  intensity  altema 
tion  with  the  progression  of  J. 

The  value  of  band  system  is  9/2,  so  the  final  expressions  for  the 

H  function  beccmie: 


16ir*L 


H 


v*v*’.  J*J**  “  9hc  Vv*\  J»J»*  J**ei 


(2.33) 


2 

where  |  R^|  q(v*v*')  replaces  to  indicate  the  lack  of  band  strength  information, 

and  Sj„j,  takes  on  the  appropriate  value  from  the  following  list: 


J*'  even,  R  branch 
J'*  odd,  R  branch 
J**  even,  P  branch 
J**  odd,  P  branch 


2(J*’  ♦  1) 

(J”  ♦  1) 

3J” 

J” 


The  Franck-Condon  factors  used  in  this  calculatioo  were  proridad  by  NiohoUs  fRsf.  21) 
and  are  baaed  on  Morse  potentials.  Thble  V  givde  these  values. 

The  value  of  |  R^  I  ^  need  in  this  calculation  was  assumed  to  be  that  correopoadii«  to 
an  f-number  of  0. 1,  solely  on  Indirect  experimental  evidence  |aee  disoussicu  by 
Gilmore  (Ref.  22)). 


Cqastioae  2.31-2.33  and  2.23  were  incorporated  into  a  digital  computer  code  which 
was  used  to  construct  a  spectral  atlas  for  this  system.  This  atlas  was  also  merged 
into  the  SACHA  air  atlas  mentioned  previously. 


it 
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Table  V,  Franck- ConckMi  Factors  For  The 

Birge-Hopfield  No 


le  No  (b»^Su  - 
•  1  ^stem 


X 


B 


B 

V”  1 

1 

2 

3 

4 

5 

6 

1 

Q 

2. 1190-7 
2. 7849-6 
1.8669-5 
8. 5095-5 
2.9664-6 
8.4338-4 
2. 0367-3 

2.9921-6 
3.4604-5 
2.0290-4 
8.0341-4 
2.4137-3 
5. 8603-3 
1.1957-2 

2.  1884-5 
2.2112-4 
1. 1238-3 
3. 8204-3 
9. 7434-3 
1.9804-2 
3.3239-2 

1.0971-4 
9. 5936-4 
4. 1742-3 
1. 1986-2 
2. 5373-2 
4. 1831-2 
5. 5152-2 

4.2151-4 
3. 1524-3 
1. 1554-2 
2. 7389-2 
4. 6537-2 
5. 9053-2 
5. 5978-5 

1.3171-3 
8.3013-3 
2. 5089-2 
4.7522-2 
6. 1416-2 
5.4340-2 
2.9782-2 

3.4716-3 

1.8101-2 

4.3828-2 

6.3144-2 

5.6360-2 

2.7371-2 

3. 0393-3 

■ 

7 

8 

9 

10 

11 

12 

13 

1 

n 

7.9072-3 
3.3321-2 
6. 2072-2 
6.  2689-2 
3.1176-2 
2.7340-3 
5. 8297-3 

1.5831-2 

5.2317-2 

7.0591-2 

4.2684-2 

5.4214-3 

5.3553-3 

2.9014-2 

2. 8210-2 
7. 0244-2 
6.2211-2 
1.4943-2 
2.1711-3 
2.8169-2 
3.4894-2 

4.5158-2 

8.0177-2 

3.8661-2 

1.4776-4 

2.2324-2 

3.8433-2 

1.3957-2 

6. 5404-2 
.  6375-2 
1.2476-2 
9. 7782-3 
4.0811-2 
2. 1321-2 
1.7964-5 

8.6174-2 

5.8193-2 

5.3382-5 

3.3398-2 

3.4640-2 

1.3800-3 

1.5049-2 

1.0372-1 

3. 1969-2 
9.5731-3 
4.6879-2 

1. 1261-2 
7.4268-3 

3. 2478-2 

■ 

14 

15 

16 

17 

18 

19 

20 

I 

n 

1. 1442-1 
8. 9762-3 
3. 3555-2 
3.6128-2 
1.1900-4 
2.9736-2 
2.2115-2 

1.1595-1 
2. 1381-6 
5.3490-2 
1.2112-2 
1.5406-2 
3.4215-2 
1.8886-3 

1.0613-1 
9. 6578-3 
5.3909-2 
2.4528-6 
3.7429-2 
1.3956-2 
7.0055-3 

9. 2899-2 
3.3955-2 
3.4531-2 
1.2953-2 
3.7755-2 
9.4965-9 
2.8850-2 

7.3580-2 
6. 2614-2 
1.0367-2 
3.8361-2 
1.5946-2 
1.4218-2 
2.9672-2 

5.3735-2 
8.4465-2 
1.0483-5 
5. 1042-2 
1.6971-4 
3. 5964-2 
7.9907-3 

3.6175-2 

9. 2570-2 
1.2100-2 
3.8781-2 
1.1410-2 
3.3118-2 
1.4963-3 

■; 

21 

22 

23 

24 

25 

KH 

27 

1 

Cl 

2. 2435-2 
8.6338-2 
4.0365-2 
1.3774-2 
3.7500-2 
9.  8604-3 
2.  2638-2 

1.2804-2 
7.0204-2 
6.9563-2 
5. 5983-5 
4.8499-2 
6.3770-4 
3.7228-2 

6.7151-3 
5.0455-2 
8.6062-2 
1. 15bO-2 
3.2247-2 
2.0551-2 
2.  1406-2 

3. 2302-3 
3.2307-2 
8.4950-2 
4.1815-2 
7. 2453-3 
4.4596-2 
1.0002-3 

1.4220-3 
1.8514-2 
7.0140-2 
7. 1932-2 
1.2946-3 
4.2193-2 
1.0148-2 

5.7133-4 
9. 5167-3 
4.9651-2 
8.5977-2 
2.3401-2 
1.6791-2 
3.8019-2 

2.0881-4 

4.3901-3 

3.0553-2 

8.0219-2 

5.8100-2 

8.9475-5 

4. 5685-2 
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2.3  THE  ©2  SCHUMANN- RUNGE  CONTINUUM 
3  -  3  — 

The  B  -  X  (Schumann- Runge)  continuum  of  O2  is  a  strongly  absorbing 

system  in  the  ultraviolet  and  is  an  important  contributor  to  the  absorption  of  heated 
air  up  to  about  10, 000*K.  Elvans  and  Schexnayder  (Ref.  23)  have  performed  a  com¬ 
bined  e>q;>erimental  and  theoretical  study  of  this  continuum  and  have  presented  results 
which  are  readily  adapted  to  our  compilation.  These  results,  which  are  listed  in 

Table  VI,  give  values  of  a  reduced  absorption  coefficient,  k ,  where  —  . 

^max 

Here  is  the  value  of  the  maximum  spectral  absorption  coefficient  of  a 

standard  sample  of  pure  O2  at  300* ,  and  K  is  the  spectral  absorption  coefficient  of 
the  same  sample  at  some  other  wavelength  and  temperature. 

Ihe  values  of  the  average  absorption  coefficient  presented  in  the  tables  for  air  in 
appendix  A  were  obtained  by  use  of  the  expression 

^SjS-Rcoot.  -  C<02>  (^XdT  •  Kav  •  "0  <*•*<> 

where  C(02)  is  the  O2  concentration  in  moles  per  mole  of  original  cold  air  and 
been  obtained  by  interpolating  linearly  with  photon  energy.  The  average 
value  of  K  was  assumed  to  be  the  value  at  the  center  of  the  energy  interval.  Since 
K  is  a  slowly  varying  ftmctlon  of  photon  energy,  the  error  introduced  may  be  neglected 
here.  The  4(>0(cm*^)  contained  In  the  above  expressloa  is  the  best  esttmate  for 
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Table  VI.  The  Abaorption  Coe£Cicient  of  Oxygen  at  Selected 
Wave  Lengths  and  Temperature 


l^oton 

Energy 

9.54 

9.18 

8. 86 

8.55 

8.43 

6.27 

(eV) 

Wave 

Length 

A  =  1.400A 

A  >  1.470A 

A«  1,500A 

TCK) 

K,  Best  estimate 

;  all  states  excited 

300 

0.000 

0.425 

0.985 

0. 995 

0.940 

0.822 

1,  000 

— 

0.424 

0.937 

0.903 

0.847 

0.744 

2.000 

— 

0.400 

0.715 

0.670 

0.590 

3,  000 

0.359 

0.674 

0.580 

0.543 

0.495 

4.  000 

0.207 

0.304 

0.559 

0.483 

0.457 

0.422 

5.000 

0.185 

0.251 

0.466 

0.404 

0.385 

0.358 

6.000 

0. 158 

0.216 

0.394 

0.343 

0.332 

0.307 

7.000 

0.130 

0. 189 

0.335 

0.296 

0. 288 

0.  270 

8.000 

0. 105 

0. 168 

0. 291 

0.260 

0.253 

0.243 

9.000 

0.081 

0. 153 

0.255 

0.224 

0.218 

10.000 

0.058 

0. 140 

0.227 

0.202 

0.200 

0. 197 

Photon 

Eneror 

8.00 

7.75 

7.47 

7.29 

7.09 

(eV) 

Wave 

Length 

X  »  1.  550A 

A  =  l.OOOA 

A  =  1, 660A 

A  -  1.700A 

A«  1,750A 

300 

0.583 

0.335 

0. 135 

0.060 

0.018 

1.000 

0.549 

0.351 

0.180 

0. 104 

0.053 

2,000 

0.471 

0.353 

0.229 

0. 170 

0.112 

3,000 

0.407 

0.330 

0.242 

0. 198 

0. 143 

4,000 

0.360 

0.  303 

0.237 

0. 197 

0. 153 

5,000 

0.316 

0.272 

0.222 

0. 189 

0. 153 

6,000 

0. 275 

0.  240 

0.200 

0. 175 

0. 148 

7,000 

0.245 

0.217 

0.180 

0. 163 

8,000 

0.221 

0.  197 

0. 170 

0.153 

0. 132 

9,000 

0.200 

0. 182 

0.158 

0. 144 

0.125 

10,000 

0.182 

0.  169 

0. 147 

0. 135 

0.119 
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Section  III 

ABSORPTION  BY  ATOMS 

3.  1  ENERGY  LEVELS  AND  OCCUPATION  NUMBERS 

In  ordei  to  calculate  photoioi’.’zation  absorption,  one  must  specify  the  set  of  atomic 
states  or  energy  levels  present  in  the  gas  under  given  equilibrium  conditions  and 
compute  the  absorption  from  each  of  these  states  individually.  The  contribution  to 
the  absorption  coefficient  from  an  individual  atomic  state  i  can  be  expressed,  for 
calculational  purposes,  as 


=  N.  y]  (T  (e)  cm  ^  (3.1) 

j  ^ 

3 

where  N.  is  the  number  density  (in  atoms/cm  )  of  atoms  in  the  state  i  ,  and  fijle) 
is  the  photoionization  cross  section. 


The  subscript  j  arises  from  the  fact  that  each  atomic  state  will,  in  general,  give 

rise  to  several  different  photoionization  transitions,  corresponding  to  different  states 

of  the  residual  atom.  For  each  of  these  different  final  states  j  there  is  a  separate 

T 

cross  section  a. .(e)  and  different  energy  threshold  e..  such  that,  to  obtain  the  total 

ijV  /  exT  jj  . 

cross  section  for  all  final  states,  the  sum  over  j  as  indicated  in  Eq.  (3. 1)  must  be 

T 

performed  for  ail  states  j  for  which  c  e.j  . 


The  theory  whereby  the  cross  sections  fjj(e)  are  obtained  is  described  in  the  next 
section.  The  present  section  is  devoted  to  a  brief  discussion  of  the  number  densities, 
or  occupation  numbers,  N^  which,  together  with  the  (T.j(e)  ,  make  up  the  absorption 
coefficient  according  to  Eq.  (3.  1). 
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The  computer  code  PIC  which  was  employed  to  calculate  the  atomic  photoiooization 
abeorpUon  coefOcienU  as  deOned  by  Eq.  (3. 1)  was  devised  to  use  mm  input  a  set 
of  occupation  number  fractions  p^  ^i^v  ^  ^  different  density  values  for 

aqy  croup  of  sUtes  i  which  one  wishes  to  consider  at  a  given  temperature  (  is 
the  total  particle  number  density.  ^  >.  For  the  code  to  calculate  the  cross  sections 

for  the  relevant  transitions,  a  spectroscopic  identification  of  the  state  i  includii^  its 
rnergy  must  be  given  with  each  set  of  6  values  of  .  Each  atomic  sUte  is.  therefore, 
spec/fied  by  an  identification  number  ii  r  n  f  8  LI  and  the  energy  of  the  state  (in  eV 
relative  to  the  ground  state  of  the  neutral  atom  of  the  species); 

•  i  is  a  two-digit  designation  of  the  species  and  charge  state:  01.  02.  etc.,  for 
oaygen  I.  oxygen  n.  etc. :  and  11.  12.  etc. .  for  nitrogen  1.  nitrogen  II.  etc. 

•  V  is  the  two-digit  designation  of  the  core,  and  is  listed  in  Table  Vn.  mkotm 
with  the  relevant  kmizatio.  notentials. 

•  n  (two  digiu)  is  the  principal  quantum  number  of  the  excited  electron  or 
outermost  populated  shell. 

a  f  is  the  orbital  angular  momentum  of  the  excited  electron  or  outermost 
populated  shell  (summation  conventions  ill  and  iv>. 

a  (28  1)  Is  the  spin  multiplicity  (summation  coni  ention  ii). 

a  L  is  the  total  orbital  angular  momentum  (summation  convention  1  and  ii). 


Table  Vn.  Table  of  Atomic  Core  Configurations  and  lonUation  Potentials* 


Species 

SUte 

Core  Designation 

o 

N 

O-I 

Is^  2s^  2p^  (^8)nl(^*L) 

V  •  1 

13.617 

(*i» 

2 

16.942 

(^P) 

3 

18.631 

2s  2p'*  (‘•p) 

4 

28.42 

5 

34. 19 

<*s> 

6 

37.88 

(^P) 

7 

39.99 

K  v> 

2p 

8 

52.  11 

*lonUatk>n  potentials  arc  given  to  lowest  sUte  of  lowest  term  of  next  stage  of  ionization. 
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Table  VII  (Cont'd. ) 


Species 

state 

Core  Designation 

jloniz.  Potentials  (eV)^| 

O 

N 

o-n.  N-I 

la^  2a^  2p^  (^P)nl(^^L) 

V  *  1 

35. 154 

14.548 

(*D) 

2 

37.659 

16.447 

(^S) 

3 

40.50 

18.601 

2e  2p^  (®S) 

4 

42.623 

20.400 

(^D) 

5 

50.028 

25.989 

(^P) 

6 

52.  79 J 

28.094 

(*D) 

7 

58.336 

32. 430 

(^8) 

8 

59.580 

33. 787 

<*P) 

9 

61.238 

35. 230 

2P*<®P) 

10 

70.339 

42.3 

(*D) 

11 

72. 127 

43.1 

<*8) 

12 

77.708 

47.4 

o-in,  N-n 

U‘*2e^2p  (*P)nK^^L) 

y  «  1 

54.946 

29.617 

2a  2p*  (^P) 

2 

63.758 

36. 703 

<*D) 

3 

70.671 

42. 145 

<*S) 

4 

75.312 

45.861 

(*P) 

5 

77.309 

47.716 

2p^(S) 

6 

83.607 

52.781 

<*D) 

7 

86.568 

54.799 

<*P> 

8 

90.642 

58.188 

o-iv.  N-ni 

le^  2a* 

y  -  1 

77.411 

47.426 

2a  2p  <’P) 

2 

87.575 

55. 757 

(*P> 

3 

97.081 

63.629 

2p*(’P) 

4 

103.881 

69. 189 

(*D) 

5 

106. 122 

70.843 

<*8) 

6 

113.088 

76.606 

O-V.  N-IV 

la*  2a  (*S)  nK*®^*L) 

y  »  1 

113.898 

77.450 

2P  <*P) 

2 

125.821 

87.442 

O-Vl.  N-V 

la*  (*S)  id(**^*L) 

y  -  i 

138. 110 

97.863 

•lonUaUoii  potentials  arc  given  to  lowest  sUte  of  lowest  term  of  next  stage  of  kmlxaUon. 
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The  Input  to  the  PIC  code  then  coneleU  of  cards  with  the  /ollowins  information: 

l.y.n.l.2S4^1.L  Energy  (eV)  P(l>.  P(2).p(3).p(4).p(5).p(6) 

The  arguments  for  p  indicate  the  different  density  values  chosen.  Since  the  number 
of  atomic  energy  levels  is  essentially  infinite,  various  conventions,  lumpings.  and 
simplifications  had  to  be  adopted  to  hold  the  number  of  levels  within  tractable  limits. 


As  the  occupation  number  calculation  depends  only  on  the  energy  of  an  atomic  state  it 
is  convenier.t  to  combine  statM  of  nearly  equal  energy  and  exploit  the  well-known  sum 
rules  over  Wigner  coefficients  in  the  evaluation  of  the  optical  transition  probabilities. 
The  following  conventions  were  adopted: 


(i)  n  >  3 

(ii)  n  2  4 

(iii)  6  £  n  <  9 

(iv)  n  >  9 

(V)  n  >  17 


L  and  S  specified  for  each  atomic  term.  Exception:  Levels 
summed  over  total  angular  momentum  L  for  given  multiplic¬ 
ity  (2S*  1)  when  experimentally  the  L-splitting  of  the  levels 
is  much  less  than  the  S-splitting  or  when  the  energies  are 
poorly  known  (Labelled:  L  »  9  ) 

Summed  over  all  L-S  terms  (Labelled:  2S  1  >  0  . 

L  »  0) 

Summed  over  electron  angular  nnomentum  1^4  (Labelled: 

I  *  4) 

Summed  over  all  values  of  the  electron  angular  momentum 
(Labelled:  I  «  0  ) 

Not  included 


The  original  calculations  performed  with  the  PIC  code  (Ref.  5)  employed  a  set  of  such 
occupation  number  fractions  obtained  from  earlier  calculations  (Ref.  24)  and  given 
(in  abbreviated  form)  in  AFSWC  TR-61-72  (Ref.  25).  Since  this  original  set  of 

numbers  covered  only  the  temperature  range  kT  »  2  cV  to  kT  »  20  eV 
it  was  inadequate  for  the  present  calculation  and  a  new  set  of  occupation  numbers 
having  the  required  format  had  to  be  computed. 
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The  thermodynemic  functions  (frscUonsl  electronic  populstions.  or  "Boltzmaw 
fractions*')  needed  to  compute  atomic  occupation  numbers,  once  the  total  atomic  con* 
centrations  are  known,  have  been  calculated  by  Gilmore  (Ref.  26).  These  calculations 
were,  of  course,  performed  for  a  particular  selection  of  atomic  enerty  levels  with 
different  identification  and  lumping  conventions  than  those  chosen  for  the  PIC  code. 


Because  of  the  existence  of  these  results  due  to  Gilmore,  our  occupation*number  cal¬ 
culations  were  confined  to  unfolding  the  Gilmore  BolUmann  fractions,  in  a  manner 
appropriate  for  PIC  input,  and  to  devisii^  a  code  (the  "EquUibrium  Code."  see 
Appendix  B).  for  computing  the  atomic  and  molecular  mole  fractions,  or  species 
concentraUons.  Thereupon,  by  mulUplylng  theee  two  reeults  (BolUmann  fracUons  x 
mole  fractions)  the  occupation  number  fractions  were  obtained  in  the  required  form. 

A  computer  code  was  devised  to  unfold  the  Gilmore  BolUmami  fractions,  attach  the 
atomic  idenUfication  required  by  PIC.  and  mulUply  by  the  speciee  concentraUons 
obtained  from  the  equilibrium  code.  This  code  caUed  "Multiplet  ID  Split"  producee. 
as  its  output,  a  set  of  cards  in  the  PIC  input  format  (previously  described)  with 
occupaUon  numbers,  fractions,  and  spectroscopic  identifications  as  requirsd. 

In  the  unfolding  process,  the  original  average  energy  of  the  combined  levels  wmm 
retained  for  each  of  the  unf<4ded  levels  and  the  occupation  numbers  were  split  in  pro¬ 
portion  to  the  statisUcal  weights  of  the  unfolded  subgroups  of  levels  of  which  the 
origins!  "level"  was  consUtuted.  For  example,  if  GUmore  had  lumped  aU  n  -  3 
levels  of  a  given  configuraUon  together,  these  were  split  by  the  fracUons 
(2f  **-  l)/2r«  (f  «  0. 1. 2)  to  give  occupation  numbers  for  the  separate  I  values. 

These  were  then  further  split  the  ratloe  (28^  *  2)/a(28^  ^  1)  and  28  /2(8  1)  to 

account  for  the  two  spin  mulUplicities  involved  when  the  older  n  -  3  electron  is 
auached  to  a  core  of  spin  8^  .  8ince  the  PIC  code  accepts  levels  sianmed  over  L- 
values.  the  Gilnnore  fracUons  were  not  split  into  subfracUons  for  the  various 
appropriate  L-values.*  As  suted  above,  the  average  energy  as  originaUy  assigned 


*8ome  of  theee  fractions  for  the  lower  levels  are  already  given  by  GUmore  for  individual 
atomic  terms,  with  L  and  8  specified  as  needed  for  PIC.  Theee  were,  of  course, 
left  as  is. 
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by  Gilmore  ik’es  rctAined  for  the  unfolded  levels  and  the  occupation  numbers  of  the 
unfolded  levels  compiled  onl*-  from  their  statistical  weights.  Thus,  the  Boltzmann 
factors  exp  (-  AE^/kT)  ,  wiaere  AEj  is  an  miergy  difference  between  the  average 
energy  E  as  assigned  by  GUmore  and  the  actual  sublevel  energy  Ej  .  have  been 
neglected.  In  addition,  this  procedure  will  cause  errors  in  the  photoelectric  edge 
positions  of  amount  AE^  .  However,  both  these  errors  should  be  less  than  the 
uncertainties  in  the  total  calculation  arising  from  other  sources.*  (It  perhaps  should 
be  emphasized  that  the  extrapolation  on  the  energies  of  the  bound  states  required  in 
the  Burgess -Seaton  theory  to  obtain  the  phase  shift  at  the  edge  is  not  carried  out  in 
the  PIC  code  by  use  of  the  input  state  energies.  It  is  done,  rather,  from  Ubular  sUte 
energies  stored  in  the  code.) 

Table  Vin  lisU  a  sample  set  of  unfolded  occupation  number  fractions  produced  by  the 
Multiple!  ID  Split  code.  The  values  listed  are  for  the  states  of  neutral  oxygen  present 
in  air  at  10.000*K  and  densities**  from  10  to  10®  normal  air  density,  denoted  by  o(l). 
a(2)  .  etc.  The  identification  index,  which  wo  have  labelled  i  in  this  section,  is 
actuaUy  carried  in  the  code  as  two  indices.  J  and  s  .  The  atomic  state  is  defined 
by  J  .  and  the  species  by  s  .  Column  1  of  Table  Vm  lisU  the  spectroscopic  identifi¬ 
cation  of  the  state.  Column  2  lisU  the  energy  of  the  stoU  (in  oV).  and  Columns  3-8 
list  the  occupation  number  fractions  Pj  »  for  the  6  density  values,  where  Nj 

is  the  number  density  of  particles  in  the  suu  i  and  N  is  the  total  number  density. 

As  indicated  above,  they  are  normalized  such  that  multiplication  of  these  occupation 
number  fractions  Pj  by  the  total  particle  number  density  ^  (electrons  - 

and  molecules,  of  course  —  are  included  in  as  separate  spcicies)  yields  the  actual 
oocupation  number  Nj  .  (A  table  of  the  values  of  for  the  required  temperatures 
and  densities  therefore  forms  part  of  the  input  information  to  the  PIC  code. )  All  the 
levels  for  which  Boltzmann  fractions  were  computed  by  Gilmore  have  not  been  included 


*  Particularly  the  matrix -element  values. 

**It  will  be  seen  that  a  second  card  has  been  inserted  after  each  primary  card  to 
carry  the  two  additional  density  columns,  since  the  PIC  code  can  accept  a  maximum 
of  six  density  columns  at  one  time. 
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Table  Vm.  Occu|Jitlon  No.  Fractlona  |)(J,i)  for  SlalM  ol  Ncultal  Otiygen  in  Air.  Col.  1  IliU  (he  Spoctioieoptc 
IdenUficalion,  Col.  2  litu  the  ener(iei,  and  CoU.  3-8  lUi  the  Occupalion  No.  FractloM 
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in  that  sUtes  of  principal  quantum  number  n  ^  5  have  been  neRlectcKl  for  some  cxcitc<l 
cores.  These  levels  have  generally  ver>'  low  population  at  the  temperatures  of  this 
study,  and  are.  therefore,  not  important. 

3.2  ATOMIC  PHOTOIONIZATION  ABSORPTIO.V 


Within  the  limits  of  the  present  ^  tudy  it  was  not  practical  to  include  the  effects  of 
atomic  lines.  The  manner  in  which  these  can  be  included  and  an  estimate  of  the  effects 
of  such  line  contributions  can  be  obtained  from  a  previous  study  (Ref.  5).  The  primar>’ 
atomic  elf^  Jt  which,  it  will  be  seen,  becomes  very  important  toward  the  high-tempera¬ 
ture  low-density  region  of  our  calculations  is  that  of  photoionization. 

A  computer  code  called  PIC  (^lotoipnization  calculation)  developed  under  the  previously 
referenced  study  (Ref.  5)  was  used  in  the  present  work  to  rompiae  the  photoionization 
contribution  (except  for  the  0“  photodetachment  contribution,  uhich  is  considered 
separately  in  Section  IV).  A  brief  summary  of  the  calculational  basis  of  this  code  (due 
to  R.  R.  Johnston,  and  taken  from  Ref.  5.  with  appropriate  modifications)  is  included 
below  for  completenees.  For  further  deUils  of  thU  calculation  and  an  actual  listing 
of  the  code,  the  reader  should  consult  Ref.  5. 

The  atomic  energy  levels  included  in  the  present  calculation  (from  which  photokmization 
****'•4^®***  t*ke  place)  and  the  concentrations,  or  occupation  numbers,  of  these  states 
have  been  described  in  the  preceding  section. 

As  stated  previously,  the  absorption  coefficient  for  photoionization  transitions  can  be 
expressed  in  terms  of  the  total  photoionization  cross  sections  <yj.(«)  from  initial 
sutes  i  to  final  sutes  J  . 
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The  general  initial  state*  (i>tklSL)  treated  in  the  present  calculation  may  be  written 


with  the  K-shell  la  understood  and  the  conficuration 


2s 


specified  1^'  the  core  index  v  (Table  VU).  A  general  photoionization  transition  from 
such  an  initial  state  may  be  written 


2»"*  Sp"**  "*(*•-)  •*  *•"*  (*'**c)  "’**  ** 


n*  n*  /S 


where  the  final  state 
momentum  r*  and  a 


of  a  fi 
idual  atom 


electron  with  energy  c  and  orbital  angular 


In  the  one-electron,  electric -dipole  approximatlor  considered  here  the  possible  final 
sutes  are  limited  by  selection  rules: 

(a)  8*  -  8  .  for  electric  multipole  transitions  in  L8  couplii^. 

(bi  L'-L-I.L.L**^!  only,  lor  dipole  transitions. 

(c)  r*  •  (^  e  I  e  0  .  where  f^  is  the  initial  orbital  angular  momentum  of  the 
^ectron  being  ejected. 


*As  defined  in  the  preceding  section,  i  and  y  are  code  numbers  representing. 
respecti%’ely.  the  ionic  species  and  core  coupling  for  the  sute  involved.  The 
re*naining  sjrmbols  n(8L  are  standard  atomic  spectroscopic  notation  for  principal 
quantum  number,  orbital  single-electron  angular  momentum,  total  spin,  and  total 
orbital  angular  momentum,  respectively. 
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For  the  photoionization  from  a  multiply -occupied  initial  aute.  a  parenUge  expansion 
(Ref.  27)  will  give  rise  in  c^neral  to  se\'eral  possible  states  of  the  residual  atom  - 
each  with  a  different  energy  threshold  and  coefficient  of  fracUonal  parentage  F 
An  enumeration  of  the  possibilities  follows:  ^ 

(a)  Ejection  of  an  outer  (n.f)  electron 

i.  •  2p  **  nl  (®l)  -  2«"*  2p"P 


Her*  n;  .  n,  .  nj,  -  Hp  .  8^  .  S"  .  Sj^  .  .  L  •  -  L,^  . 

f*  *1^1  2T  0  . 

(b)  Ejection  of  a  2p  outer  electron:  Decompoeing  the  initial  state  into  parenU: 


1 


every  pair 


for  which 


K 


p  •'SL 


n  -I  \  2 

P  <8pLp> 


is  nom'anishing  constitutes  a  possible  final  state  of  the  residual  atom  with 
"i  ‘  "s  •  "p  •  "p  “  ^  •  Sc  •  ®p  •  ^c  "  S  -  0  or  2  . 

<c>  Ejection  of  a  2s  outer  electron:  Decomposing  the  Initial  state  into  parents. 
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the  only  possible  states  are 

n  =2:  S=0.L*0;S^=^,L=L*0.F  *2 
•  P  ^  P  p 


or 


"s  S=2*^*®**p*®**-p**-*®*Fp=l. 

in  both  of  which 

"1  *  **.  “  ^  •  *  L-.  .  *•  =  1  • 

s  s  c  p  c  p 

(d)  Ejection  of  a  2p  inner  electron.  Decouplinc  one  2p  electron  from  the  rest 

2.  •  2p  P  ^  ‘*^12)  (**-)  -  2  *P  [*•"*  ^  (*%)  *P  (****-12)  "*  (**•)] 


■  IP 


8,„L, 


X2*'l2  11*^  ^®P^P^  **®12^12 


)! 


the  outer  (n.l)  electron  must  be  coupled  to  /  ^^p)  * 
initial  state  may  be  written 


so  the  recoupled 


Where  -  FJ^U^(IL^Ll  ;  L’Xjg)  U*(|  8^8  |  ;  S^Sj^)  and  the  U-huKrUons 
are  the  Jahn  coetaclenta  of  Ref.  28.  Experimentally,  the  enersies  of  sutes 
with  different  values  of  L*')  usually  vary  much  less  than  do  the  enerkiea 

of  sutes  with  different  ^  ^^p)  *  Thmrwtorm  thm  (^*L")-spllUins  of  the 

photokmization  edges  is  ignored  and  the  sum  over  (^*L**)  is  carried  out 

n  n  -I  /8  \ 

directly,  leading  to  final  sUtes  2s  2p  ^  f  *’l  1  nf  with  probability 

F*  .  So  n*  ■  n.  .  n*  -  n„  -  I  ,  8^  •  S^  .  L  -  .  T  -  O.  2  . 

P  sspp  epep 
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(c)  Ejection  of  a  2s  inner  electron.  Decomposing  the  initial  state  as  in  Case  (d) 


F  (S' 


(I  V  2  =  '*••*12)  •  •'p  - 


As  in  Case  (d).  the  sum  over  S**  is  carried  out  directly,  leading  to  final 
states  2s  •  2p  P  (***1-12)  prohabUity  .  So  n^  «  n^  -  1  . 

”p  *  "p  *  ®c  *  ®p  •  ^c  *  ^12  *  F  -  I  . 


Such  a  decomposition  of  the  initial  state  into  substates  each  with  a  different  ionization 
potential  resulu  in  a  splitUng  of  the  photoionizatkm  edges.  For  calculations,  o.g. . 
of  Planck  mean  opacities  at  high  temperatures, this  splittU^  should  yield  resulu  only 
slightly  different  from  a  calculaUon  in  which  spliUing  U  ignored.  However,  the  split¬ 
ting  of  the  2p-edgee  is  tjrpicaily  a  few  electron  volu  in  magnitude  and  therefore  should 
affect  the  present  low -temperature  absorption  coetficlenu  significantly. 


Oure  the  initial  atomic  sute  i  Is  decomposed  into  all  possible  final  sUtes  of  the 
residual  ion.  the  photoionization  cross  section  is  evaluated  ui  the  one-electron, 
electric -dipole  approximation.  For  a  transition  from  an  iniUal  bound  electronic 
state  (n.f)  to  a  final  state  of  the  free  electron  with  enorg>'  <  .  the  photolonizatlon 
cross  section  may  be  written  (Ref.  24) 


3  *> 

.  .  I6r  e**^ 
*  3c 


(3.3) 


for  photon  lrequenc>’  .  In  terms  of  the  integrals  o\'er  the  radial  ^^‘ave 

functions 


t 

I 


/  r"  dr  |,4r)  rR^^r)  (3.4) 

This  result  has  been  derived  sssuminc  L*S  couplinc  with  the  neglect  of  spin-orbit 
intemcUon.  and  assumins  that  the  wave  function  of  the  final  free  electron  state. 

R^,(r)  .  is  normalized  on  the  enercy  scale 

/  r^  dr  R^,|,(r)  R^.(r)  -  -  c)  (3.5) 

Thus  the  only  unknown  quantities  in  Cq.  (3. 3)  are  the  redial  matrix  elements  R|^  . 

In  most  previous  work  (Refs.  24.  29.  and  30)  these  matrix  elements  have  been 
approximated  by  hydrocenlc  values. 

For  electrons  with  initial  state  principal  quantum  number  n  2  6  or  orbital  ancular 
momentum  1^2  hlfdrocenic  reeults  have  been  used  here.  For  smaller  valuee  of 
these  quantum  numbers,  some  recently  obtained,  and  significantly  more  accurate 
expressions  for  the  radial  matrix  elements  are  used.  For  small  values  of  the  elec¬ 
tron  final  eneny.  the  approximation  of  Burgess  and  Seeton  (Ref.  32)  is  used  and  is 
discussed  in  the  first  of  the  subsequent  subsections.  An  approximatioo  valid  for 
high-energy  ^ectrons  is  described  in  the  second,  and  the  third  presents  the  method 
whereby  the  two  approximations  were  joined. 

3.  3  THE  LOW-E.N'ERGY  THEORY 

Recently,  burgess  and  Smiton  (Ref.  32)  have  presented  an  approximation  to  the  radial 
matrix  cimnents  in  terms  of  the  asymptotically  correct  wave  functions.  This  approxi¬ 
mation  derives  from  the  observation  of  Bates  and  Oamgaard  (Ref.  33)  that  the  major 
contributiofi  to  the  radial  integral  for  bound-bound  transitions  usually  comes  from 
values  of  r  sufficiently  large  that  the  effective  potential  is  a  Coulomb  potential. 
Replacing  the  actual  one-electron  wa\’e  function  by  its  asymptotic  form  -  a  linear 
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combination  of  the  regular  and  irregular  Coulomb  wave  function  for  the  ubser\’ed 
value  of  the  energy,  modified  for  amall  r  to  ensure  convergence  of  the  radial 
integrals  -  Bates  and  Damgaard  evaluated  the  radial  matrix  elemenu  and 

presented  their  results  in  tabular  form. 

®*^C*^*  Seaton  <ltef.  32)  applied  similar  com ide rations  to  the  e\'aluation  of  the 
radial  matrix  elements  lor  bound-free  transitions.  Wliereas  the  asymptotic 

behavior  of  the  bound -state  wave  function  U  determined  the  physically  obscr\'ed 
energy  of  the  bound  state,  the  large-radius  behavior  of  the  free-electron  wave  function, 
at  a  given  energy,  is  determined  a  phase  shift  c)  .  In  the  approximation 

of  aaymptoUcally  correct  wave  functions.  Burgess  and  Seaton  numerically  e\'aluatcd 
the  radial  matrix  elements  and  parameterized  the  resulting  photoionization  cross 
section  |Eq.  (3.3))  in  the  form 


tic) 


V  ^  l^l.r.*^;c)  ♦  A,*(L*,c)J 


t*\ 


1 1  ♦  cv*) 


•*V||.(c)-l 


(3.6) 


ahero  L*  is  the  orbital  angular  momentum  of  the  toul  system  in  the  final  suto.  The 
quantities  C|^,(L')  are  the  reeulU  of  the  angular  integraUons  and  may  be  simply 
expressed  in  terms  of  Racah  coefficienU  aa  described  in  Ref.  32 .  G  .  y  .  and  ^ 
arc  Ubulated  by  Burgees  and  Seaton  aa  functions  of  for  various  combinations  of 
(  and  r  .  and  <  is  the  electron  kinetic  energy  in  rydber«is  divided  by  .  The 
basic  variable  of  the  theory*,  v  .  is  defined  by 


^  (lA"  -  ia) 

where 

1/A^j  -  1  By  cm"^  -  109.737.3  cm"* 

l/X  -  term  value  of  the  initial  state  (cm"S 


40 


1/V 


t«rm  limit  obtained  removing  the  initial  (n.f>  electron  and 
leavinc  the  residual  ion  in  iu  final  state  (cm'S 
charce  of  the  residual  ion 


The  dintensional  constant  a  of  Eq.  (3.6)  U  given  in  terms  of  the  Bohr  radius  a  In* 


4r  2  /  V*  \ 

»  oiz^) 


O*-) 


(3.8) 


where  e(v)  is  a  normalization  correction  to  the  asymptotic  bound-sUte  «rave  fiMKrtion 
and  may  be  estimated  from  knowledge  of  the  physical  bound  states  of  the  system.  For 
example,  the  ev^uatlon  of  t(v)  lor  the  outer  electron  in  the  configuration 

2s  *  2p  ^  j  nl  is  obtained  from  consideration  of  the  series  obtained  by 

varying  n  with  everything  else  fixed.  The  observed  energies,  r  ...  of  the  terms 

n  I 

in  the  resulting  series.  relaUve  to  the  series  ionization  limit  1.  -  Urn  c  .  are 
written  "  • 


<Vl  -  V  ■  *  <*  •> 

in  terms  of  their  effective  quantum  number  .  and  the  quantum  defect,  .  is 

defined  by 


n*f 


(3.  10) 


Seaton  (Ref.  34)  has  shown  that  the  so-deliaed  quantum  defect  p  .  is  a  continuous. 

m  - 

anal3rtic  function  of  the  energy  (for  fixed  f  )  which  may  be  written  p|(-  l/v*)  . 
Further.  P|('  l/v^  must  tend  to  an  integer  as  e  tends  to  any  one  of  the  Integers 
0  .  I . f  .  Then  riri  nf  Cq.  (3.0)  is  given  by 


-  I 


•  I 


di* 


(3.  II) 


e\'aluftted  at  the  observed  physical  enero*  of  the  state  tnf)  whose  normalisation 
is  required. 

Thus,  to  estimate  the  obser\’ed  energy  of  at  least  one  additional  tern  of  the 

(n*0  eeriM  of  the  initial  sute  (n.O  is  required.  For  outer  (n.  f)  electrons,  such 
sutes  are  usually  known  -  at  least  upon  relaxation  of  the  condition  i'-.r  the  outer 
^^^^oplin®  L  be  fixed.  For  inner  electrons.  howe\*cr.  in  the  presence  of  an  excited 
outer  electron  [e.R.  .  the  2p  orbital  in  the  state  2s^  2p*  (^p)  3s  of  .Nl)  .  the  excited 
sute  necessary-  to  estimate  Im  a  doubly  excited  sUtc  -  2s"  2p(“p)  3s  3p  in  the 

preceding  example.  The  energies  of  such  states  are  generally  not  known.  The  present 
calculation  uses  an  approximation  of  Seaton  <Ilef.  34)  for  such  inner  2p-orbitals 

«  lil-"  -  2) 

and  chooses  c  *  1.0  for  inner  2s «oibitals.  Some  comparisons  with  Hart ree-Fock 
solutions  indicate  that  this  approximation  is  reasonable. 

For  the  phase  shift  t)  of  the  free-electron  final  state.  Seaton  (Ref.  34)  has 

shown  that  a  low-energy  approximation  is  gi%’en  by  cxtrapolatiiqt  the  quantum  defect 
P|(<)  .  Eq«  (3.  10).  of  the  series  containing  the  final  state  to  small  positive  emorgies 
according  to  the  formula  (cf.  also  Ref.  35) 

£|,(«.L*)  *  x’P|«<€)  (3.12) 

For  electrons  with  principal  quantum  number  n  >  3  .  the  initial -sUte  oceufiation 
numtiers  arc  summed  owr  the  total  orbital  angular  momentum  of  the  atom,  the 
decomposition  Into  final  1.*  values  is  therefore  ignored,  and  the  phase  shift  is  taken 
to  tie  zero  and  f^/(2f  •  1).  where  U  the  greater  of  f  and  f  .  Since 

eUi'trons  with  orlilul  angular  momentum  f  ^  2  in  the  final  state  have  small  iihasc 
shifts  at  low  energies,  the  presemt  calculation  again  assumes  the  phase  shift  to  vanish. 
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For  photonionlzation  of  jui  ouler  electron  to  a.  final  s  or  p  state  the  zero  etncrty 
please  shift  is  obtained  Eq.  <3. 12).  extrapolating  linearly  from  the  physical  boimd 
states  with  high  principal  quantum  numbers.*  Howe\’er.  in  the  photoionization  of  an 
inner  electron  in  the  presence  of  an  older  excited  electron  a  difficulty  occurs  similar 
to  that  encountered  for  such  states  in  the  normalization  procedure.  Again,  the  physi> 
cal  states  required  are  doubly  excited  states  and  are  unknown;  the  extrapolation  of  the 
quantum  defect  is  therefore  impossible.  The  present  calculation  sets  the  phase  shift 
to  zero  for  such  states.  Generally,  the  most  significant  low-energy  phase  shift  is 
for  final  s-states.  thus  the  2p  cs  photoionization  cross  section  in  the  presence  of 
an  outer  excited  electron  is  in  question.  However,  the  2p  -•  cd  transition,  with 
small  d-state  phase  shift,  usually  dominates  the  2p  —  cs  .  so  this  uncertainty  is 
pnfoably  not  too  great.  A  more  serious  difficulty  is  the  neglect  of  the  p-state  phase 
sluft  in  the  2s  —  cp  transition  in  the  presence  of  the  excited  outer  electron.  The 
uncertainty  due  to  this  approximation  can  probably  only  be  resolved  by  detailed  electron 
scattering  calculations  or  Hartree-Fock  calculations  of  the  continuum  state. 

A  low-energy  expansion  of  the  Burgess -Seaton  cross  section.  Cq.  (3.6).  for  integer 
V  and  to  first  order  in  c  agress  with  a  corresponding  expansion  of  the  exact  hydro- 

genic  results.  Thus,  the  Burgess  -Seaton  theory  is  to  be  considered  a  valid  approxi- 

2 

matlon  when  the  electron  kinetic  energy  is  much  less  than  Z  rydbergs.  Furthermore, 
when  sufficient  information  regarding  the  physical  bound  stales  is  available,  the  result 
of  Burgess  and  Seaton  is  probably  more  reliable  than  continuum  Hartree-Fock  calcu¬ 
lations  (Refs.  36  and  37)  as  some  effects  due  to  exchange  and  polarization  of  the  core 
are  reflected  in  the  physical  energies  of  the  states  with  large  principal  quantum 
number  used  in  the  extrapolation  for  the  phase  shift. 

Thus  we  conclude  that  the  Burgess -Seaton  approximation  is  a  reasonably  valid  and 
useful  approximation  near  the  photoionization  threshold,  except  for  the  aforementioned 
inner-shell  transitions  in  the  presence  of  an  excited  outer  electron.  The  simplicity  of 


*The  bound-state  energies  required  are  tabulated  in  the  computer  code  so  that  high 
accuracy  is  not  required  of  the  input  state  energies. 
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Ui©  resultii^  expression  lor  the  photoionixation  cross  section  renders  it  particularly 
suitable  for  e\*aluation  of  the  larse  number  of  initial  states  present  in  a  calculation  of 
this  t3rpe. 

3.4  THE  HIGH >E.\ERCY  THEORY 

For  energies  of  the  final  state  free  electron  much  jnneeter  than  zero,  the  Buriecss- 
Seaton  approximation  is  inapplicable.  An  approximation  of  frequent  utility  in  hish- 
enercy  acatterin®  calculations  is  the  Bom  approximation.  Hou'evcr.  a  straiithtlonibarti 
applicaUon  of  the  Bora  approximation  to  calculation  of  the  photoionUalion  cross  sec¬ 
tion  leads  to  an  incorrect  result,  as  discussed  in  Ref.  5. 

The  basic  expression  for  a  radiative  transition  -  in  the  electric  dipole  approximaUon  - 
betareen  atomic  states  |i>  and  |f>  is 

2  2 

"  ”~c  *  I<^l'r|i>I*  <3.13) 


from  which  Eq.  <3. 3)  is  obuined  upon  performing  the  ansular  Integrations  and  the 
average  over  the  polarUation  direcUons  of  the  incident  .'adiation.  As  shown  by  Bethe 
and  aalpeter  <Ref.  3t).  forms  equivalent  to  Eq.  <3. 13)  are  given  by 


|<flPH>l* 

m  cw 


4t^V 


<3.  M) 


These  are  the  dipole -velocity  and  dipole-acceleration  forms,  respectively.  Z  is  the 
nuclear  charge,  and  m  the  electron  mass.  It  can  bo  shown  (Ref.  5)  that  a  use  of  the 
acceleration  form  of  the  matrix  element  in  a  high  energy  Born  approximation  gains 
one  iterate  of  the  Bora  series  o%'er  the  use  of  the  dipole-^-elocity  form.  Accordingly, 
the  present  calculation  e%’aluales  the  photoionization  cross  section  for  high-clectron 
crerg>*  b>*  the  dipole-acceleration  Born  approximation. 
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When  the  angular  intecrationa  and  polarization  averacea  are  performed  in  Eq.  (3. 14), 
a  reault  analocoua  to  Eq.  (3.3)  ia  obtained,  namely. 


aj^<c)  = 


16r^Z^e^ 

3m^u;^c 


(3.  15) 


with 


<3.  16) 


In  the  first  Bom  approximaUon  the  properly  normalUed  free  ■  electron  wave  f^tkm 
is  the  plane  wave  -  (k/2rS^^^  «*P  (OT  •  7)  -  so  that  Eq.  (3. 16)  becomes 


'W'*  <*• 

o 

in  terms  of  the  spherical  Beesel  functions  j^kr)  .  where  -  electron  kinetic 
enency  in  rydbergs. 

Through  uee  of  a  generalization  of  the  HF8  code  of  Hermann  and  fliillman  (lief.  39). 
wave  functions  have  been  generated  for  the  many  atomic  atatee  preeent  in  the  gas. 
The  reeulting  numerical  wave  funcUona  have  been  fitted  by  analytic  functions  of  the 
form 


iViin 


2 


-o_r 


(3.  m 


45 


with  intac^r  valuer  of  .  From  comparison  of  the  numerical  results,  the  foUowinfs 
conventions  were  adopted  in  order  to  reduce  the  number  of  orbitals  required  while 
stai  retaining  accuracy  suitable  to  the  present  application: 

(a)  For  ground>state  cores  -  2s  *  2p  P  nf  .  the  orbital  (n.O  was  calculated 
forSsnsSandf^  0. 1.2.3  only.  For  all  values  of  (n.l)  the  core 
orbitals  are  approximated  by  those  calculatcsd  with  (n  <  3.1)  unless  1  =  3. 
when  the  core  orbitals  used  are  those  calculated  with  (n  *  4  ,  f  «  3)  . 

(L)  For  excited  cores  with  one  or  more  2s-electrons  excited  to  a  2p-orbit.  the 
orbital  (n.l)  is  calculated  for  n  »  3.4.5  and  I  «  0.1,3  and  for  n  =  3.4 
for  1*2.  For  larger  n  the  (n.l)  orbital  is  approximated  by  those  cal¬ 
culated  in  the  presence  of  the  ground-state  core.  The  core  orbitalc  arc 
subject  to  the  same  conventions  as  in  Case  (a). 


In  terms  of  the  bound  orbitaU.  Eq.  (3. 18).  Eq.  <3. 15)  for  the  cross  section  becomes 


(3.  19) 


where 


and 


m 

l|(o./l)  -  f  s’®*  *^J|<*)  d* 


r(. .  §) 


(1  oS 


l/2(f'»^'»l) 


(3.  20) 
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These  are  most  easily  evaluated  uainc  the  recursioo  relations 


l|4a.^)  *  ^0*11  ”  ^  ♦  I  -  1)  -  1) 

-  IVMl  ♦  sin  13  tsn‘*  (l/o)l 

1^0. 0)  -  -  o  *)  l,.|<o.O) 

lg<a.O)  •  tsn**  (i/or) 


(3.21) 


As  discuaaed  by  Chandrasekhar  (Ref.  40)  the  three  expreaalona  |  Eqa.  (3. 13)  and 
(3. 14)1  are  equivalent  for  exact  wave  functions  but  for  approximate  wave  (unctions  the 
fact  that  they  weight  different  regions  of  configuration  space  differently  leads  to 
different  results.  Thus  the  Burgees  -aeeton  approximation  of  wave  functions  correct 
for  large  r  starts  from  Cq.  (3. 13)  and  conforms  to  the  famllia4  jncertainty  principle 
arguments  that  low  ensrg)f  scatteruig  states  preferentially  sample  the  l<^-range 
parts  of  the  bound  state.  flimUarly  the  high-energy  scattering  state  weights  more 
heavily  the  short-range  part  (the  hlgh-momentum  Fourier  components)  of  the  bound 
state,  as  does  the  acceleration  form  of  the  matrix  elemen*..  Variatioiially  determined 
arave  fieictioas  are  usually  less  reliable  for  small  than  for  intermediate  values  of  r  . 
so  the  accuracy  at  short  distances  of  the  HF8  wave  funcUona  represents  a  source  of 
uncertainty  in  the  present  results. 


3.S  THE  APPRQXIMATC  PHOTOIONIZATfON  CIIO0B  SECTION 


The  occurrence  of  the  spherical  Bessel  function  in  Eq.  (3. 17)  reprseents  an  approxi¬ 
mation  to  the  more  correct  positive  energy  regular  Coulomb  function  valid  if  k  »  Z  - 
a  domain  of  validity  opposite  to  thid  of  the  Burgess -Soeton  approximation,  as  dis¬ 
cussed  at  the  end  of  the  first  subsection.  For  the  present  calculation,  an  exponential 
interpolation  formula  is  used. 

etc)  •  ■  0*^*“)  (3.21*) 
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^B8  ^HE  photoioni ««t  ion  croM  s^cttaos  cmlculittad  accontti^  to 

tkm  Biirgc— »8— ton  and  high-oneny  appronimationa .  raapactivaly.  Tha  paramoCar  a 
araa  clwaaa  to  waishi  tiM  two  appttaOmatioaa  agnally  at  c  »  rydbarga. 

Compariaona  of  tkm  plntolonlaidion  croaa  aacUon  lor  OI  and  NH  TTtrtalafid  from 
Cq.  (3. 22)  witk  Hartraa'-Fock  calcalationa  cbtalnad  aritli  tlm  compatar  program  of 
Oalgarao.  Hanry.  and  fltawart  (Raf.  26)  and  witli  Iqrdrogonic  raanlu.  ara  ptaaotad  la 
M-  2  and  arill  not  ba  rapaatad  bara.  It  la  aaaa  from  thaaa  compariaona  that  tba 
bydroganic  approaimaUon  lor  tba  2p  adgaa  conalatantly  undaraatlmataa  tba  Hartraa- 
Fbck  raaalta  and  agraamant  improraa  for  atataa  of  highar  dagraa  of  kmixatlon.  TIraa. 
tba  bydroganlc  appcoaimatloo  la  moat  aultabla  at  higb  tamparaturea  and  low  dmialtiaa. 
Comparing  continuum  Planck  maana  obtalnad  by  fltawart  and  I^att  (Raf.  29)  with  tboaa 
obtalnad  In  Rid.  S  baara  out  this  obaarvatlon. 

To  account  for  bound-bound  tranalUona  with  uppar  level  above  the  merging  limit,  tba 
pbototonlsation  croaa  aectlon  la  extrapolated  linearly  to  a  value  (zVn^)  rydbarga 
below  the  vacuum  edge,  wdiera  n^^  la  the  principal  quantum  number  vdiare  merging 
of  the  llnaa  occura.  Thia  la  dlacuaaed  further  In  Ref.  5  where  the  formula  for  n^ 
la  given. 
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8«ctl<»  IV 

MISCELLANEOUS  CONTRlBUnONS  TO  THE  ABSORPTION 


4. 1  NO2  •  ABSORPTION 

4.1.1  Exp^rlmMital  Dftta 

Mueller  (Ref.  •)  hue  eompUed  the  available  experlmeatal  data  oa  the  abaorptloa  coef- 
flcleat  of  NOij*  Siaoe  experlmeatal  data  exiet  oaly  ia  a  very  narrow  raace  of  wave* 
leastba.  1000*K  as  T  ^  2000*K.  Mueller  extrapolated  the  data  to  cover  the  rai^ 

2500  A  s  X  s  5000  A  and  300*K  s  T  s  S000*K.  Ualng  the  raeaaurementa  of  Wurater 
and  Marrone  (Ref.  41).  QIbnore  (Ref.  42)  corrected  aome  of  theae  predicted  valuea  and 
gave  a  prediction  up  to  8000*K. 

In  our  application  we  uae  the  valuea  of  Mueller  aa  corrected  and  enlarged  by  Gilmore. 
(See  Figure  1  and  Table  DC. ) 


4.1.2  Calculation 

The  absorption  coefficient  M^02  ^  ^  ^  calculated  within  the  range  of 

parameters: 


1000*K  s  T  s  8,000*K 
Range  I  10  a:  p/p^  ^  10"® 

19,  837  A  s:  X  2:  1,167  A  (in  linear  steps  of  0. 1  eV) 

where  the  temperature  and  the  relative  density  p/p^  of  the  air  sample  are  provided 
by  the  equilibrium  code  (Appendix  B). 


49 


7000 


Table  DC.  NO^,  •  AbaorpUoe  CoefDciMit  to  cm  /% 


Wave- 

length 

X 

Cn) 

¥ 

1.9837 

5041 

1.4188 

7058 

1.1020 

9074 

.5834 

.5220 

.4723 

.4312 


.3873 

.3420 


0.3006 
0.2834 
0. 2681 
0.2543 
0.2419 


2 


23191 
2 


33267 

35286 


3 

3 


41339 


1000 


R 


1.4l 

7.5^ 

1.4* 

6.0* 

1.7* 

3.4* 

4.9* 

6.2* 

6.7* 

6.4* 

5.0* 

4.0* 

3.3* 

2.8* 

2.5* 

2.2* 

o  n3 


0 

0 

0 

2.  ft 
2.3* 
3.6* 
1.1* 
2.2* 
3.2* 
4.0* 
4.6* 
4.9* 
4.6* 
4.1* 
3.7* 
3.4* 
3.2* 
3.0* 
2.83* 
2.6* 


O 

0 

l.ll 

1.4* 

4.2* 

6.0* 

1.5* 

2.3 

2.8* 


3.38* 

3.40* 

3.35* 

3.30* 

3.25* 

3.20* 

3.16* 

3.12* 

3.10* 


O 

2.8* 

1.6* 

4.9* 

8.2* 

1.0* 

1.7* 

2.25* 

2.7* 

3.0* 

3.15* 

3.25* 

3.2* 


3.0* 

2.94* 

2.85* 

2.80* 

2.75* 

2.70* 


2.4* 

2.4* 

2.4* 
3 


2. 

2.4* 
o  ^3 


0 

0 

0 

0 

2.0* 
2.0* 
2.0* 
2.0* 
2.0* 
2.0* 
2.0* 
2.0* 
2.0* 
2.0* 
2.0* 
2.0* 
2.0* 
2.0* 
2.0* 
2. 


The  experimental  values  of  the  absorption  coefficicmt  k  are  given  in  cm  /gm  of  NOn 

-1 

From  this  we  can  find  the  absorption  coefficient  In  cm  by  multiplyii^  by  the 

mass  of  NO2  per  cm^: 


#^N02  =  •  ^<N02)  •  M(N02)  •  ^ 


(4.1) 


5 


witfc  the  aKtlMiilar  waIgM  M(N(^)  -  M.Ol  g  mole**,  the  gas  comtam  R  -  S2.0S3  atm 
•*•>'*  Btole**  aad  tka  temperatare  T .  The  laole  tnctkm  X(NO^)  aad  the  iveMare  P 
emm  b*  fir—  from  tte  •qnlllbrliim  ood*. 


TMm  DC.  tiM  experlmeaCal  daf  for  «  cover  tlie  rmi«e  of  paranieCers: 

1000*K  s  T  2  8000*K 

R—fe  n 

1R.8S7  A  X  ^  S41R  A 

la  that  part  of  Itaas*  I*  which  to  aot  cowered  fay  Raage  O.  to  set  to  aero. 

To  ftod  the  ral—e  of  «  for  raluM  of  X  faetwe—  the  X-grto  of  Table  DC.  the  follow- 
lag  laterpolattoa  nethod  to  appUed.  Let  aad  X,  be  two  tabulated  valu—  of  X  la 
Table  DC,  aiHaoeat  to  a  glirea  X  •  Thea  r  {X)  to  fouad  by  uae  of  the  expr— aioa 


»«I*<C(X)  -  logK<Xi)  ♦  (logictXg)  -  log«(Xij^^^*  <4.2) 

4.2  O'  PHOTO-DETACHMENT 

4. 2. 1  Experimeatal  aad  Theoretical  Data 

The  crou  aectioo  for  0“  photodetachmeat  up  to  3  eV  haa  beea  noeaaured  by 

Smith  (Ref.  43).  For  thia  eaergy  raage  only  the  tranaitioa  to  the  ^P-oaygen  level 
with  the  threshold  energy  1.456  eV  contributes.  The  energy  region  above  3  eV  .  whereg 
in  addition,  transitions  to  the  and  levels  contribute  with  threshold  energies 
of  3.43  eV  and  5.66  eV  respectively,  has  not  yet  been  investigated  experimentally. 
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TlMoreticml  derlvmUoM  of  tke  croos  sectkai  <?PD  ^  nuide  by  Klebi  aad  Bruedaer 
(Ref.  44)  aad  Martto  aad  (hooper  (Ref.  4S)  do  aot  0ve  results  m’lUdi  mstcli  tlie  expert^ 
mental  vahiss  in  tks  energy  region  above  S  eV.  The  values  calculated  in  Ref.  45  are 
assumed  to  be  accurate  only  within  §0~~70  pmrcent. 

Smith  (Ref.  43)  gave  a  prediction  of  the  crons  section  SpD  for  the  range  3  eV  ^  he  £ 
7. 5  eV.  He  assumed  that  the  cootrSiutlons  of  the  ^P.  ^D.  ^8  components  are  pro* 
port  tonal  to  the  product  of  statistical  weight  and  photon  energy  according  to  Bates  and 
Massey  (Ref.  44).  In  addition,  be  used  the  sscperlmcntal  threabold  energies. 

For  our  application,  we  took  the  data  measured  (s  3  eV)  aad  predicted  (e  3  eV)  by 
Smith.  For  7.5  eV  n  he  ss  10.425  eV.  we  esctrapolated  Smith's  curve.  (See  Figure  2 
and  Table  X.) 

4.2.2  Calculation 

The  absorption  coefficient 


^PD  “  "PD  ■  **v(0') 


(4.3) 


is  calculated  within  the  range 

1000*K  £  T  £  24,000*K 
10  a  p/p^  a  10"® 

19.837  A  2  X  2  1.167  A 

(0.625  eV  <  hi/  :s  10.625  eV  in  linear  steps  of  0. 1  eV) 

where  the  temperature  T  and  the  relative  density  p/p^  are  provided  by  the  equilibrium 
code. 
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Figure  2  Photodetachment  Cross  Section  of  O 
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TiUde  X.  CroM  SectSoa  for  O*  rHotodetmchmet 
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N^(0-)  «  Hy  .  XCO-)  (4.4) 

U  m  pmkiet  of  the  totml  Mtmber  daeelty  Ny  akSid  tl»  mole  fractkm  X(0“)  .  geaeraied 
by  the  eguiUbrliim  code. 

Ummmr  taterpotethm  Is  applied  to  flml  for  the  vahaee  of  hr  .  aot  given  to 

Table  X.  This  yields  the  expressloa 


"PdO")  “  W-5) 

where  hr|  and  hr2  are  the  Uhular  points  of  Table  X.  adjacent  to  hr  . 

4.3  FREE -FREE  ABSORPTION 
4.3. 1  Basic  Equations 

We  consider  an  air  sample  consisting  of  neutrals,  ions,  and  electrons.  Let  NyCZ)  and 
NyCe)  be  the  number  density  of  Z-fold  charged  positive  Ions  and  of  electrons,  respec¬ 
tively  for  a  given  sample  in  equilibrium,  these  number  densities  depend  on  the  tem¬ 
perature  T  and  the  relative  density  p/p^  and  may  be  taken  from  the  equilibrium 
code.  The  positive  ions  are  assumed  to  have  a  pure  Coulomb  field  of  charge  Z  . 

The  absorption  coefficient  Ppp  in  cm"^  for  free-free  transitions  of  the  electrons  in 
the  field  of  the  positive  ions  can  be  written  as  (Refs.  7,  30): 

^FF  I  Nv(e)Ny(Z)gK(Z)  •  (4.6) 

Z 
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Krmxnwr  a  semSclaMicml  absorfsUoa 
em^  in  Eq.  (4. 6>  Is  gjLvmn  bgr 


co«IEIcl«il  KiX}  jpmr  ninctron  pnr  ctn  pnr  Um 


«»  ■  dw  {^)  #  (# 


1/2 


<4.T) 


mmd  g  is  thn  Gsmt  fsolor  av«nc«d 
Boltemsnn  vslocllgr  dfattrlbiitios  of  tiM 
is  cn*^.  ssd  m  Um  msss  of  tkm 


ths  (soo^rslstivistic. 

TIm  qsstiHqr  V  is  tiM 
tiM  ociMr  plusrsiesl 


Tbs  GsuBt  fiictor  i  as  s  fsnetioii  of  ths  dimsasioslsss  gusaUUss 


is 


_ iv 

kT 


2  2 

r  -  z*  • 


■  (R^:  Kjrdbsrg  constant) 


(4.8) 


<4.») 


has  boon  givsn  by  Kanas  and  Lattsr  (Rsf.  47). 


4.3.2  Calculation 


-X 


Ths  absorption  coofUciont  Mpp  bi  calculatsd  for  ths  range  of  parameters 


Range  I 


1,000* 

10 


s  T 


p/p. 


24.000*K 

io-« 


19, 837  A  2  X  *  ^  2  1, 167  A  (in  linear  stepa  of  0. 1  eV) 
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where  the  temperature  T  and  the  relative  density  p/p^  are  provided  by  the  equilib¬ 
rium  code.  The  necessary  physical  constants  and  conversion  factors  are  (Chemical 
Rubber  Handbook.  1962—63): 


The  number  densities  Ny(e)  ,  Ny(Z)  can  be  fouivl  from  the  total  number  density 
and  the  mole  fractions  generated  by  the  equilibrium  code  by  use  of  the  formulas: 

Ny(e)  =  •  K(e)  (4. 10) 

N^l)  =  Ny  •  ^(^2)'^  X(CO'^)  +  XCNO*^)  +  XCN"^)  +  X(0^)  +  XCC"*^) 

*  X(A*)  +  X(Ne*)|  (4.11) 

Ny<2)  -  Ny  •  +  X<0**)  +  XiC**)  *  XiA**)  +  X(Ne**)J  (4. 12) 

where  the  sums  in  (4. 11)  and  <4. 12)  are  taken  over  all  single  and  double  charged  ions. 

A  two-dimensional  array  of  Gaunt  factors  g  (u^V^  over  the  range 

-  3.  10.  100,  1000 

Range  n 

0.25  £  u  s  30 
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is  taken  from  Karzas  and  Latter  (see  Table  XI).  The  u-steps  are  chosen  such  that  the 
following  linear  interpolation  of  g  as  a  function  of  log  u  and  log  y  can  be  made. 

Let  .  (“a^2)  ’  (^*^1)  ’  ^  points  of  the  (ur^-grid  of 

Table  XI  which  are  adjacent  to  an  arbitrary  point  (uy^)  of  Range  n  (Figure  3). 


The  first  interpolation  step  results  in  the  formula 


*K)  -  +  |«(“b>'f)  -(«  Vf)l 


log  -  log  Uj 


with  i  «  1,  2  .  From  these  values  »  the  final  value  gi(uY^  is  obtained: 


(4. 13) 


g(iir^  -  +  |,«(wr2)  - 


2  2 
log  y  -  log 

log  yg  -  log  yj 


(4. 14) 


In  that  part  of  Range  I,  which  is  not  covered  by  Range  n,  /Xpp  is  set  to  zero. 


li 
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Table  XI.  Gaunt  Factors 


=  3 

o 

II 

=  100 

y^  =  1000 

0.  2 

1.50 

1.38 

1.  20 

1.  100 

0.3 

1.425 

1.325 

1.  16 

1.090 

0.4 

1.375 

1.285 

1.  14 

1.070 

0.6 

1.320 

1.240 

1.  13 

1.060 

0.8 

1.280 

1.23 

1.  11 

1.055 

1.0 

1.265 

1.20 

1.  105 

1.050 

1.2 

1.25 

1. 1925 

1.  100 

1.048 

1.4 

1.235 

1. 1825 

1.095 

1.047 

1.6 

1.225 

1. 1770 

1.0925 

1.045 

1.8 

1.215 

1. 170 

1.090 

1.044 

2.0 

1.2075 

1. 166 

1.086 

1.0425 

2.4 

1. 195 

1.  160 

1.085 

1.0413 

2.8  . 

1.185 

1.155 

1.082 

1.040 

3.2 

1. 176 

1. 150 

1.080 

1.0390 

3.6 

1.  170 

1. 1475 

1.079 

1.0380 

4.2 

1. 160 

1. 1450 

1.0785 

1.0370 

5.0 

1.150 

1. 140 

1.079 

1.0365 

5.8 

1.140 

1. 1375 

1.0795 

1.0375 

6.2 

1.135 

1. 1350 

1.0798 

1.0380 

7.2 

1. 125 

1. 1340 

1.080 

1. 0390 

8.2 

1. 115 

1. 1325 

1.0805 

1.0395 

9.2 

1. 105 

1. 1300 

1.081 

1.0400 

10.2 

1.095 

1. 128 

1.0815 

1.0410 

12.0 

1.078 

1.125 

1.084 

1.0420 

14.2 

1.060 

1. 120 

1.0855 

1.0435 

16.2 

1.045 

1.117 

1.0870 

1.0450 

18.0 

1.030 

1.114 

1.0885 

1.0460 

20.6 

1.015 

1. 1075 

1.0905 

1.0465 

23.0 

1.000 

1. 1020 

1.0910 

1.0475 

25.4 

0.985 

1.0975 

1.0925 

1.0485 

26.4 

0.980 

1.0950 

1.0936 

1.0490 

28.0 

0.970 

1.090 

1.0940 

1.0495 

30.0 

0.960 

1.086 

1.0950 

1.0500 
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Section  V 
DISCUSSION 


5. 1  GENERAL  REMARKS 

A  survey  of  the  tables  leads  to  the  observation  that  the  region  of  the  temperature- 
density  plane  covered  by  these  calculations  cannot  be  readily  partitioned  into  subregions 
where  either  atomic  transitions  or  molecular  transitions  are  clearly  dominant.  In 
fact,  our  results  indicate  the  existence  of  a  lai^e  region  where  both  atomic  and 
molecular  transitions  contribute  appreciably  to  the  absorption,  with  the  relative  con¬ 
tribution  depending  on  photon  energy  (or  the  equivalent  frequency).  Certain  general 
conclusions,  however,  can  be  drawn  on  the  basis  of  the  calculation: 

<  1)  Atomic  photoionization  contiibutlons  are  more  important  at  lower 
temperatures  than  has  been  freqpiently  assumed 

(2)  Ultraviolet  molecular  bands  of  N2  are  very  Important  up  to  about 
IS.OOiTK  at  the  higher  densities 

(3)  Much  more  experimental  work  and  analyses  are  needed  In  all  areas 
covered  by  the  calculation. 

5.2  COMPARISON  WITH  MEYEROTT.  80K0L0FF.  AND  NICHOLL8 

Considering  the  differences  between  the  present  work  and  that  of  MSN  (Ref.  1).  the 
agreement  between  them  is  certainly  acceptable.  As  an  example,  a  graph  from  Ref.  1 
Is  presented  In  Figure  4  with  our  results  for  the  total  air  absorption  plotted  thereon. 

In  our  resultSySome  of  the  prominent  band  characteristics  exhibited  by  their  curve 
are  suppressed,  with  the  largest  disagreement  centered  on  the  N^  1*  at  about 

3. 1  eV.  Our  curve  rises  above  theirs  both  In  the  ultraviolet  and  in  the  Infrared. 

The  former  effect  Is  due  to  enhancement  of  the  ultraviolet  absorption  by  the  Blrge- 
Hopfield  bands  of  N2  which  were  not  Included  in  the  MSN  compUatlon.  The  difference 
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Figure  4  Absorption  Coefficient  of  Heated  Air  vs.  Photon  Energy  (T  =  12,000*K.  p/p  =  1) 
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in  absorption  in  the  infrared  is  due  to  the  larger  nitrogen  photoionization  given  by  the 
present  computations. 

5.3  THE  ATOMIC  PHOTOIONIZATION  CONTRIBUTION 

# 

A  very  brief  comparison  with  experiment  of  the  PIC  code  results  was  performed  in 
Ref.  5.  In  the  context  of  the  present  work,  it  is  instructive  to  compare  also  with  the 
previous  tabulation  of  MSN  and  with  the  more  recent  calculation  due  to  Ashley  (Ref.  48). 

The  photoelectric  contribution  was  neglected  in  the  MSN  calculation  for  temperatures 
T  <  12,000*’K.  Since  12,000*’K  is,  in  fact,  the  upper  limit  of  their  tabulation,  com¬ 
parison  is  possible  only  at  this  one  temperature.  Figure  5  shows  this  comparison  for 
nitrogen.  The  centers  of  the  squares  represent  the  MSN  points,  and  the  centers  of 
the  circles  show  the  results  of  the  present  calculation  utilizing  the  PIC  code.  It  should 
be  noted  that  the  PIC  code  result  lumps  all  contributing  stages  of  ionization,  whereas 
the  MSN  calculation  included  only  the  contribution  of  the  neutral  atom. 

The  simple  hydrogenic  model  utilized  by  MSN  is  seen  to  overestimate  the  present 
results  (except  at  the  very  shortest  wavelengths)  by  an  amount  varying  from  agree¬ 
ment  to  about  a  factor  of  four.  This  is  somewhat  reminiscent  of  the  comparison  of  the 
unit -gaunt  factor-hydrogenic  aiid  I^C  results  in  Ref.  5,  where  the  hydrogenic  result 
for  the  2p-edges  was  also  shown  to  be  high  at  low  i^oton  energies  (~  4  eV  or  less) . 

In  tie  present  case,  since  we  are  considering  more  edges,  the  behavior  is  more 
complicated,  but  toward  X  ~  0. 14  p  and  less  the  PIC  results  have  increased  over 
the  MSN  hydrogenic  values.  Just  prior  to  the  lower  wavelength  limit  of  the  calcula¬ 
tion,  the  lowest  nitrogen  2p-edge  enters  the  range,  which  causes  the  absorption 
coefficient  to  shoot  up  maxicedly. 

Also  shown  on  the  figure  are  the  results  of  Ashley  for  the  more  limited  wavelength 
range  X  =  0.38  to  0.65  p.  The  remaining  figures  show  further  comparisons  with 

-4  • 

Ashley's  results.  Figure  6  is  for  nitrogen  at  6,000*K  and  relative  density  p/p^  =  10 

Figures  7  and  8  show  oxygen  at  T  =  8,000*’K  and  at  relative  densities  p/p.  =  1 

-6  ° 
and  10  ,  respectively.  Finally.  Figure  9  compares  the  present  results  with  those 

! 
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WAVELENGTH  (MICRONS) 

Figure  8  Linear  Spectral  Absorption  Coefficient  of  Air  vs.  Wavelength 
(T  =  8000*K  .  p/p^  =  10-®) 
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of  Ashley  at  24.000”K  and  relative  density  10  .  These  conditions  have  been  chosen 

to  span,  to  a  limited  extent,  the  range  of  conditions  covered.  It  should  be  noted  that 
Ashley  lists  separately  the  contribution  from  each  stage  of  ionization,  rather  than 
lumping  them  together  as  in  the  PIC  results.  Ashley's  calculations,  like  tliose  of 
MSN,  were  based  on  the  hydrogenic  cross  sections  with  unit  gaunt  factor.  However, 
Ashley  employed  a  strong  simplification  over  the  ordinary  hydrogenic  approximation 
in  that  he  treats  the  principal  quantum  number  n  as  a  continuous  variable.  His 
photoelectric  cross  section  is  therefore  evaluated  as  though  it  aiose  from  a  continuous 
distribution  of  photoelectric  edge  strength,  and  is  correct  (hydrogenically  speaking) 
at  the  thresholds,  or  edges,  but  not  in  between.  His  cross  section  has  the  behavior 
shown  in  Figure  10  where  the  correct  hydrogenic  edges  are  indicated  and  Ashley's 


Figure  10  Photoelectric  Cross  SsciUm  (Sketch) 


cross  sscUon.  ths  dotted  line.  Interpolates  between  them.  The  fllllng-in  between  the 
edges  has  ths  effect  of  including  a  "line  contribution**  and  is  at  least  partly  desirable. 
However,  this  additional  contribution  enters  in  an  uncontrolled  fashion.  The  extra 
oscillator  strength  thereby  introduced  is  not  correlated  with  the  correct  acfclitional 
oscillator  strength  needed,  and  is.  in  fact,  divergent.  (One  can  vary  this  by  iKHing 
that  /<^v)  dv  would  be  infinite  for  \shley's  crosa  section.)  On  the  other  hand,  this 
defect  is  partly  compensated  by  the  fact  that  Ashley  concomitantly  trsaU  the  occigia- 
tkHi  numbers  as  continuously  distributed,  but  with  the  proper  totals.  This  has  the 


effect  of  reducing  the  absorption  coefficient  Na  in  partial  compensation  for  the 
leased  oscillator  strength  which  has  been  introduced.  The  general  behavior 
indicated  in  our  diagram  for  o(v)  .  however,  still  persists  in  p(i/)  even  after  this 
compensation,  as  can  be  seen  in  Figures  5  through  9.  One  would  expect  Ashley's 
method  to  be  good  for  large  n  where  the  edges  are  closely  spaced.  However,  the 
spectral  regicxi  (X  =  0. 38  to  0. 65  p)  of  his  calculation  is  primarily  the  region  between 
the  n  =  2  and  n  =  4  edges  of  N  and  O,  and  these  edges  are  quite  widely  spaced. 

5.4  RECOMMENDED  IMPROVEMENTS  FOR  FUTURE  COMPILATIONS 

Some  of  the  less  prominent  bands  of  the  Schumann-Runge  and  the  Ng  first  negative 
systems  were  not  included  in  the  calculations  because  of  the  lack  of  band  strength 
information  for  those  bands.  There  were,  consequently,  too  few  lines  for  these  sys¬ 
tems  contained  in  the  compuUtion  (see  Table  n  ).  Omission  of  these  lines  has  the 
effect  of  producing  fictional  "windows"  in  the  associated  spectrum.  In  the  case  of 
the  Og  Schumann-Runge  system,  this  leads  to  the  lack  of  a  tabulated  average  absorpUon 
coefficient  from  5. 1  to  5.5  eV.  For  Og  in  air  this  gap  is,  of  course,  filled  in  by 
•b®  jfption  arising  from  other  mechanis.ns.  F*uture  compilations,  however,  should 
include  the  weaker  bands. 

In  addition,  there  is  a  growing  concern  over  possible  inaccuracies  in  the  band  strengths 
due  to  the  use  of  Franck -Condon  factors  based  on  the  Morse  potential.  Franck-Condon 
factors  based  on  more  realistic  potentials  are  now  provided  with  the  aid  of  larger 
digital  computer  programs  and  so  those  factors  should  be  utilized  (Refs.  49  and  50). 

Also,more  of  the  ultraviolet  band  systems  of  Ng  should  be  included  in  future  compila¬ 
tions,  particularly  when  better  spectroscopic  constants  and  band  strengths  become 
available.  Gilmore  (Ref.  22)  attributes  most  of  the  Ng  absorption  in  this  region  to 
the  Birge-Hopfield  No.  2  system.  His  argument  is  based  on  an  experiment  by  Wray 
and  Teare  (Ref.  51).  We  feel  that  both  Birge-Hopfield  systems  should  be  included 
and  that  their  relative  importance,  as  well  as  that  of  other  Ng  systems,  needs  further 
investigation.  It  should  be  noted  here,  however,  that  in  problems  of  practical  interest 
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connected  with  heated  air  absorption, the  actual  ultraviolet  absorbing  mechanism  is  of 
interest  only  as  a  matter  of  detail ,  since  any  band  system  covering  the  same  energy 
interval  as  that  of  either  Birge-Hopfield  system  causes  the  air  to  be  relatively  opaque 
in  that  interval. 

In  the  case  of  the  photoionization  calculation,  the  most  significant  improvement  that 
might  ultimately  be  desired  would  be  to  assign  energies  of  spectroscopic  accuracy  to 
the  input  energy  levels.  As  noted  in  Sec.  3. 1,  the  average  energies  used  by  Gilmore 
have  been  employed,  and  this  approximation  limits  the  spectral  accuracy  of  the 
results.  For  further  discussion  of  the  limitations  of  the  theory  involved,  the  reader 
should  consult  Ref.  5. 
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Appendix  A 

AN  EQUILIBRIUM  CODE 


A.  1  INTRODUCTION 


An  important  step  in  the  calculation  of  the  radiative  properties  of  high-temperature 
air  is  a  determination  of  the  particle  concentration  of  those  species  which  contribute 
significantly  to  the  radiation  processes.  In  general,  at  the  temperatures  under  con¬ 
sideration,  air  can  be  a  mixture  of  molecules,  atoms,  and  ions  which  are  formed  by 
various  chemical  reactions  from  the  low -temperature  composition  normally  comprising 
the  terrestrial  atmosphere.  By  adopting  thecriteron  that  the  composition  of  heated 
air  can  be  fully  determined  by  the  methods  of  chemical  equilibrium,  it  is  possible  to 
establish  the  required  specie  concentrations  in  terms  of  two  independent  thermo- 
d3mamic  properties. 

Beginning  with  Brinkley  (Refs.  A-1  and  A-2),  several  methods  have  been  developed 
for  calculating  the  equilibrium  composition  of  a  chemically  reacting  system .  Of 
these,  Brinkley  adopted  the  technique  of  defining  his  system  in  terms  of  the  independ¬ 
ent  parameters  temperature  T  and  pressure  P  .  On  applying  the  principles  of 
chemical  thermodynamics,  he  resolved  his  solution  into  a  set  of  algebraic  equations 
which  on  linearization  were  solved  using  an  iteration  procedure.  Using  this  same 
method  as  a  basic  technique,  Fickett  (Ref.  A -3)  developed  a  computer  code  that  later 
evolved  as  part  of  a  more  general  HUG  code  (Ref.  A-4)  which  has  been  used  exten¬ 
sively  for  solving  shock  and  detonation  wave  problems.  Hilsenrath,  Klein,  Woolley 
and  Sumida  (Refs.  A -5  and  A-6)  used  a  different  set  of  independent  properties  in 
calculating  the  equilibrium  composition  and  thermodynamic  properties  of  a  gaseous 
system.  Selecting  the  parameters  temperature  T  and  density  p  ,  a  solution  was 
obtained  without  linearization  by  adopting  a  method  based  on  the  convergence  scheme 
of  Wegstein  (Ref.  A  -7). 


77 


In  essence,  the  determination  of  the  composition  of  a  reacting  gaseous  system  varies 
as  to  the  selection  of  independent  properties  and  the  method  used  to  solve  the  resultant 
set  of  nonlinear  simultaneous  equations.  For  the  present  study,  the  composition  of 
high -temperature  air  has  been  defined  in  terms  of  temperature  and  density.  The 
solution  of  the  algebraic  relations  has  been  iteratively  solved  using  the  techniques 
similar  to  those  of  Refs.  A-6  and  A-7,  and  in  addition  several  organizational  features 
of  the  HUG  code  have  been  adopted.  The  derivation  of  the  equilibrium  relations  in 
terms  of  T  and  p  will  be  described  using  a  generalized  notation.  In  turn,  the 
method  of  solution  will  be  indicated  and  the  particular  application  for  the  composition 
of  high -temperature  air  will  be  discussed. 

A.  2  DESCRIPTION  OF  THE  CHEMICAL  SYSTEM 

Let  the  system  be  defined  in  terms  of  s  ditferent  gaseous  species  Y^^^  ,  which  may 
be  formed  from  c  elements  according  to  the  chemical  reaction: 


c 

k  =  l 


tk 


x(*^) 


1  ^  t  <  s 


The  chemical  formula  of  Y^^^ 


is  written  symbolically  as: 


(A.l) 


"tl  "t2 


tk 


‘^tc 


1  ^  t  <  s 


(A.  2) 


Once  the  set  of  elements  X^  ^  is  selected,  applicable  to  a  particular  gaseous  system, 
the  chemical  formula  of  the  specie  Y^  can  then  be  described  by  the  vector 


(0^1  ,0^2  • 


a 


tk 


(A.  3) 


The  components  of  all  y^  vectors  form  an  o  matrix  which  provides  a  means  of 
defining  the  elemental  composition  of  all  species  Y^^^  . 
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It  is  now  possible  to  choose  a  set  of  c  ^  c  independent  vectors  [corresponding 
to  the  set  of  independent  species  and  assume  an  order  of  the  such  that 

the  first  c  species  are  the  independent  ones.  In  this  manner,  the  remaining 
dependent  sftecies  can  be  written  in  terms  of  the  independent  species  according 

to  the  reaction: 


Y^^^  =  ^  c  +  1  <  i  s  s  (A. 4) 

j  =  l 

Expressing  the  species  Y^^^  and  Y^^^  in  the  above  reaction  by  the  reaction  (A.  1), 
we  find  a  relation  between  the  a  and  v  matrices  as 


c  1  <  k  £  c 

"ik  "  X  *'ij"ik 
«  IJ  JK  c  +  1  s  i  s  s 

It  is  noted  that  only  for  the  case  c  -  c  are  the  equations  ( A.  5)  independent.  If 
c  >  c  ,  it  is  necessary  to  omit  the  extraneous  equations.  In  this  case,  it  is  possible 
to  introduce  c  independent  elements  and  (c  -  c)  dependent  elements. 


Let  n^  be  defined  as  ^he  mole  number  of  the  t-th  specie  at  some  prescribed 
thermodynamic  state  defined  by  T  and  p  .  In  addition,  let  n  be  the  total  number 
of  moles  of  gas  at  T  and  p  .  The  chemical  composition  of  the  system  can  then  be 
described  in  terms  of  the  mole  fraction  or  the  mole  number  per  original  moles 
Ct 


and 


where 


where 


(A.  6a) 


(A .  6b) 
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Here  is  defined  as  the  total  number  of  moles  at  selei^ted  standard  conditions  of 
and  .  Further  reference  to  and  c^  wUl  more  generaUy  denote  the  con¬ 
centration  of  the  t-th  specie. 

The  size  of  the  chemical  system  can  be  specified  either  in  terms  of  the  total  amounU 

of  elements  or  in  terms  of  the  set  of  independent  species.  Thus,  for  the  elements 

x(k) 

we  have: 


Qj^  GRAM  ATOMS  OF  THE  ELEMENT  X^*^^  1  <  Ic  <  c 

SirnUarly.  using  the  independmt  species  the  size  of  the  system  can  be  expressed 

as: 

Qj  MOLES  OF  INDEPENDENT  SPECIES  1  s  J  s  c 

where  the  system  is  considered  to  be  composed  of  only  the  independent  species.  The 
quantities  and  are  thus  seen  to  be  related  by  the  expression: 

Qfc  “  1  =  k  s  c  (A.7) 

It  is  noted  that  and  q^  are  conserved  in  the  course  of  all  possible  chemical 
reactions  that  may  occur  in  the  system .  For  a  given  set  of  Q|^  .  we  can  obtain  the 
number  of  moles  n  at  T  and  p  by 


"o  =  J 


(A.  8) 


The  set  is  determined  from  the  known  elemental  or  molecular  composition  of 
the  gaseous  mixture. 
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A.  3  THERMODYNAMIC  RELATIONS 


As  a  basic  assumption  of  the  analysis,  we  assume  that  all  reaction  products 
exhibit  the  characteristics  of  an  ideal  gas.  Thus,  the  ideal  gas  law 


PV  -  nRT  (A.  9) 

holds  and  relates  the  pressure  P  ,  and  the  volume  V  to  the  absolute  temperatur  .' 

T  .  R  is  the  molar  gas  constant  and  n  is  the  number  of  moles. 

The  average  molecular  weight  M  of  the  gaseous  system  can  be  defined  using  the 
molecular  weights  M^  of  the  sr>ecies  as 

s  s 

“  '  n  2  "t”t  =  f 

t=l  t  =  l 


where  the  composition  in  terms  of  c^  has  been  expressed  using  Eq.  (A.  6b). 

From  the  basic  definition  of  mass  density  p  and  the  above  relation  for  M  ,  the 
volume  of  the  system  V  is  written  as 


V  _  nM  _  %  ^ 

~  P  ~  P  t  ."l 


Thus,  Eq.  (A.  9)  may  be  rewritten  in  terms  of  p 


P  = 

M 


c^Mj 


giving 


(A.  11) 


(A.  12) 
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For  convenience,  the  compressibility  factor  Z  is  often  introduced  into  the  equation 
of  state.  This  quantity  is  defined  by  the  relation 


Z 


(A.  13) 


where  is  the  molecular  weight  of  the  unreacted  gas  mixture  at  staitdard  conditions. 
Now  the  mass  of  a  closed  system  is  conserved  when  chemical  reactions  occur.  This 
is  stated  by  the  expression 


nM  =  n  M 
o  o 

Hence,  the  compressibility  factor  is  also  given  as 


(A.  14) 


Z 


(A.  15) 


and  denotes  the  increase  in  the  number  of  moles  in  the  system  due  to  the  occurrence 
of  dissociation  and  ionization.  The  equation  of  state  may  therefore  bo  alternately 
expressed  in  the  following  form 


p  .  i^ZRT 


(A.  16) 


The  particle  or  number  density  is  often  used  to  express  the  specie  concent  ations  of 
gaseous  mUtures  and  has  the  units  of  particles  per  cm^.  Using  Eqs.  <A.9)  and 
<A.  ll)  and  introducing  Avogsdro's  number  N®  .  the  total  particle  density  of  the 
system  N  }u  given  by 


PNf  n^ 

RT  ’  V  *  M 


<A.  17) 
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In  a  similar  manner,  using  the  definition  of  the  mole  fraction  ,  the  specie  particle 
density  N  ^  is  found  to  be 


=  XjN 

(A.  18) 

A  more  convenient  expression  for  the  total  particle  density  is  found  by  taking  the 
ratio  of  this  quantity  at  T  and  p  to  that  at  standard  conditions.  Use  of  Eqs.  (A.  13) 
and  (A .  17)  gives  the  expression 

N(To.p„)  ■ 

(A.  19) 

N(Tq,P^  ia  reoagnlzad  as  Losehraidt's  number  L  .  Thus.  Eq. 

(A.  19)  becomes 

N  •  N(T.p)  -  Zl-^(p/p„) 

(A.  20) 

From  Eq.  (A.  18),  the  specie  number  density  in  terms  of  the  concentratloM  X|  and 
C|  Is 

(A.  21a) 

(A.  21b) 

where  the  relation  c^  -  ZX^  la  found  from  Eqs.  (A. 6a)  and  <A.6b). 

A.  4  EQUILIBRIUM  COMPOSmON 

For  a  gaseous  system  in  chemical  equilibrium,  the  concentration  of  the  dependent 
species  and  c^  can  be  obtained  from  the  Independent  speclee  X^  and  Cj  by 
APfdi^lnf  the  low  of  mass  action  to  the  reactions  as  given  hj  Eq.  (A. 4). 


S3 


=  “l  /7  <A.22a) 

J  =  1 

c  *►  1  ^  i  s  s 
c 

®i  °  77  <A-  22b) 

J  =  » 

The  equilibrium  consUuiU  k*  and  kl*  are  defined  in  terms  of  the  standard  free 
energies  for  each  component  in  the  reaction.  These  quantities,  expressed  in 
alternative  form  suitable  for  either  of  Eqs.  <A.  22).  nre: 

*"  “*  ‘  2  pj  5]  -  S  *  "u)  -  »]  *»  P  (A.  23.) 

*"  ‘  [•'M  RtJ  -  RT  *  [(,?,  •'o)  -  ‘J  *“  (^)  2*b) 

The  mass  balance  equatkms.  based  on  the  conservation  of  q^  moles  of  Ind^wndent 
components  .  are  written  as: 


"j  *  1^1  Vt  ’  *•)  *  »  j  s  c  <A.24) 

Oi  introducing  the  concentrations  X|  and  C|  of  Eqs.  (A. 6).  the  mass  balance 
€X|uations  become 


S4 


I 


X, .  j:  .„x, 


'j* 


£ ,  V 


I  *C^l 


X  a  J  a  c 


i 


(A.  25a) 


(A.  25b) 


The  companion  relationa  Eqa.  <A .  6)  can  alao  be  reexpreaaed  aa 


Z  X,  -  » 

t  »l  * 


(A.  26a) 


t  -i  * 


(A.  26b) 


The  solution  to  the  problem  of  determlninc  ^  equilibrium  composition  of  a  reactii« 
gaa  mixture  for  the  formulation  presented  here  is  most  readily  accorapllshsd  by 
adopting  the  relations  la  terms  of  the  concentration  c^  .  It  ie  noted  that  the  depend¬ 
ent  specie  concentrations  ehoem  in  Cq.  (A.  2Sb)  can  be  elimiaeted  by  Eqe.  (A.  22b). 
Thus,  the  problem  is  reduced  to  finding  the  nonnegative  set  oC  independent  specie 
concentretionB  c^  from  the  set  of  c  roaas  belaaoe  relatione  given  by  Sq.  (A.  25b). 
Once  the  c^  are  knosm.  the  c^  are  readily  determined  and  the  related  component 
number  deneities  and  oompreeaibUity  factor  are  than  caicniable. 

The  standard  free  energies  1^  are  generally  available  in  tabular  form  as  found  in 
the  literature.  Fbr  practical  problems.  It  is  sufficiently  accurate  to  represent  these 
data  hy  an  empirical  relation  mm  a  fancUoe  of  tempereture.  The  representation  of 
used  in  this  work  is  given  by 

C|  D|  A  C|  - 

^  -  Ilf-  ♦  Aj(l  -  to  T)  -  BjT  - 
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^1®  represents  the  standard  heat  of  formaUon.  ^hich  is  the  increment  in  enthalpy 
associated  ^ith  the  reaction  of  forming  the  given  compound  from  its  elements,  with  each 
substance  in  its  thermod>'namic  standard  state  at  the  given  teinperature.  The  coeffi- 

cienU  Aj  ,  ,  c^  ,  ,  and  are  determined  by  the  cur\'e  fitting  of  tabulated 

data. 

A.  3  EQUILIBRIUM  CODE 

The  application  of  the  equilibrium  relations  previously  discussed  are  laciliuted  by 
certain  organizational  features.  The  code  offers  provisions  for  investigating  different 
chemical  subsystems  which  are  parts  of  a  larger  chemical  system.  The  overall  sys¬ 
tem  can  be  considered  to  be  composed  of  r  different  species  which  include 

the  above  deacribod  s  spaclee  and  aa  a  sidiaet.  Thta  aubeet  can  then  be 

defined  by  ni^  and  ni|  •  which  are  the  m  niunbera  of  the  aid  i^  componenta. 
With  this  nomenclature,  the  following  identities  can  be  stated: 

*  Y  1  s  J  s  c  (A.  28) 

Hi 

r  ■  Y  c  ^  1  s  i  s  s  (A.  2D) 

Similarly .  a  subset  of  c  elements  can  be  associated  with  a  larger  set  of  o 

olemenu  with  the  identity 


1  s.  k  s  c 


<A.  30) 


\%ith  these  identities,  the  description  of  the  chemical  system  given  in  Section  B.  2 
holds  for  e%'cry  subset  as  well  as  for  the  total  system.  Thus,  all  the  relations  of 
that  section  are  %'alid  for  the  quantities  characte4*ized  by  s  bar  shown  above.  In 
addition,  the  related  quantities  fu  ..sh  the  further  identities 


o 


ik 


8C 


(A.  31) 


1 


%  =  (A.  32) 

(A. 33) 

The  use  of  this  orgsnizstkmsl  festure  in  the  equilibrium  code  is  obtsined  by  providing 
for  two  sets  of  input  dsts.  The  first  set  comprises  the  date  which  define  the  total 
overall  a3rstem  of  the  r  species  and  e  elements.  Individual  cases  are  then  specified 
by  additional  sets  of  input,  each  of  which  defines  the  thermodynamic  state  by  the 
independent  variablee  T  and  p  and  the  subeeCs  of  gaseous  components  and  elements 
of  the  subsystem. 

The  iteration  scheme  used  to  determine  the  c  independent  specie  concentrations  c^ 
has  been  adapted  in  part  from  the  work  of  Refs.  A-6  and  A-7.  It  has  been  fotsad  that 
convergence  is  facilitated  by  selecting  the  Xj  *s  which  iviU  have  dominant  concen¬ 
trations.  Thus,  in  individual  cases,  various  subsystems  can  be  defined  at  desired 
temperatures  and  densities  which  eaclude  species  that  may  not  exist  or  for  which 
thermodynamic  dau  are  incomplete.  In  turn,  this  allows  convenient  changes  in  the 
selected  Independent  species  during  aqy  one  computer  run. 

The  high-temperature  composition  of  air  determined  by  the  equilibrium  code  for  this 
study,  is  based  on  the  molecular  composition  of  dry  air  compiled  by  Gluekauf  and  as 
given  by  Gilmore  in  Ref.  A-a.  These  data  are  shown  in  Table  A-I  and  have  been 
used  to  calculate  the  for  the  total  system. 

Table  A-1.  Molecular  Composition  of  Dry  Air 


Component 

Mole  Fercent 

N 

78. 084 

O 

20.848 

A 

0.834 

C 

0.033 

Ne 

0.003 
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The  thermodynamic  data  used  to  define  the  standard  free  energy  of  the  gaseous 
components  have  been  obtained  from  several  sources.  The  values  of  the  heats  of 
formation  used  and  their  data  sources  are  listed  in  Table  A-II.  Additionally, 

the  data  sources  are  also  gdven  from  which  the  coefficients  shown  in  Eq.  (A .  27)  have 

been  derived.  The  34  gaseous  species  shown  in  this  table  comprise  the  full  set  of  r 
species  . 
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Table  A -n.  Thermodynamic  Data 


Component 

Heat  of  Formation. 
Al«,t  »‘0*K 
(kcu/mole) 

Reference 
(Appendix  A. 

Sec.  A.  6) 

Reference 
(Appendix  A, 

Sto.  a.  6) 

1.  COg 

-93.965 

9 

8.  9 

2.  NO2 

8.734 

9 

8.  9 

3.  O3 

34. 739 

9 

8.  9 

4.  CO 

-27. 190 

9 

8.  9 

5.  CO+ 

295.990 

8 

8 

6.  Ng 

0.0 

9 

8.  9 

7.  N$ 

359.310 

8 

8 

8.  NO 

21.480 

9 

8.  9 

9.  NO'*^ 

234. 880 

8 

8 

10.  O2 

-10.400 

8 

8 

11.  O2 

0.0 

9 

8.  9 

12.  Oi 

277.900 

8 

8 

13.  £• 

0.0 

8 

8 

14.  C 

169.990 

8 

8,  9.  11 

15.  C*^ 

429.832 

8 

8.  9.  11 

16. 

992. 125 

8 

8.  9.  11 

17.  C**^ 

2096.240 

8 

8.  10.  11 

18.  N 

112.535 

8 

8,  9.  11 

19.  N^ 

448.051 

8 

8.  10.  11 

20.  S*^ 

1130.970 

8 

8.  10.  11 

21. 

2224.980 

8 

8.  10.  11 

22.  O' 

25.310 

8 

8 

23.  0 

58.989 

8 

8.  9.  11 

24.  O**^ 

373.037 

8 

8.  10.  11 

25.  0*^ 

1183.770 

8 

8.  10.  11 

26.  0^ 

2450.970 

8 

8.  10.  11 

27.  A 

0.0 

8 

8.  11 

28.  A"*^ 

363.436 

8 

8.  10.  11 

29.  A*^ 

1000.530 

8 

8.  10.  11 

30.  A*^ 

1943.970 

8 

8.  10.  11 

31.  Ne 

0.0 

8 

8.  11 

32.  Ne*^ 

497.311 

8 

8.  10.  11 

33.  Ne*^ 

1444.720 

8 

8.  10.  11 

34.  Ne*^ 

2914.790 

8 

8.  10.  11 

0f 


I 
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A.  7  DIGITAL  COMPUTER  PROGRAM  LISTINGS 

The  following  list  is  the  Fortran  compile  listing  and  input  data  for  the  equilibrium 
code. 
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EQUILIBRIUM  CODE 

DIMENSION  YBC 10) .CUEB( 10) • BETA! 10) *MS(50 ) fYA(50) tALF( lOtSO) 
DIMENSION  ETH(4*50) *DTHC4*50) *CTHf4.50) *BTH(4*50) * A TH ( 4 • 50 ) *EK TH ( 4 
lt50) •MZ2(4*50) 

DIMENSION  FZ (50) vENUC 10*50) tEKC  50) *CUEC 10 ) *C ( 50 ) * CB ( 10 ) vAJJPC 10*10 
1)*A(10*10)*B(10) tCPISO) *MC(50) *OA(9*54) * TS ( 50 ) *X ( 50 ) *LS ( 10 ) *MB(50) 
DIMENSION  DEL(50) •ITER(9) *KEY( 50) •  C0NVG(50)*  CBCP(50)*  SUMT(50) 
DIMENSION  PCC(50)»  0P(9*54) 

DIMENSION  CPLAST(IO)*  CBLAST(IO)*  QA(IO)*  00(10) 

EPS*l«0E-06 
EPS2*l.E-20 
OCONVG  «  100.0 
LIMIT  *  300 
REWIND  15 
C  INPUT  I 

READ  INPUT  TAPE  5*2213 

1  READ  INPUT  TAPE  5 *  1000 *NES *NSP* ( YB ( I ) *  I « 1 *NES ) 

3  READ  INPUT  TAPE5 . 1 900 * SSWl * SSW2 * SSW3 • SSW4 • SSW5 
I-NSP4.4 

READ  INPUT  TAPE  5* 1001  * ( OA( 1 • J ) • J«5 *  I ) 

DO  650  J  «  5*  I 
650  0P( 1 • J)  »  OA( 1 tJ) 

READ  INPUT  TAPE5*1010* (CUEe( I ) *1*1 *NES) 

READ  INPUT  TAPE  5 • 1010 * ( BETA ( I ) *  I «1 *NES ) 

DO  10  I«1*NSP 

10  READ  INPUT  TAPE5* 1020*MS( 1 ) *YA( I ) * ( ALF( L *  I ) *L«1 *NES ) 

READ  INPUT  TAPE5*1030*R*ENL*RHOZ*TBZ 
DO  5  1*1*4 

READ  INPUT  TAPE  5* 1 100 * ( ETH ( I  * J ) •DTH( I  * J ) * CTh ( I  * J ) * BTh (  I  * J ) * ATM (  I  * 
IJ) tEKTHI I *J) vHZZl I tJ) tJ*! *NSP) 

5  CONTINUE 
C  INPUT  II 

100  CONTINUE 

READ  INPUT  TAPE  5* 1200 • IJX • JKX vT vOT *RHORZ *ORH 
1F( IJX-400 1801 *801 *800 

800  WRITE  OUTPUT  TAPE  15*2211*  I J:; 

END  file  15 

RFWTND  15 

CALL  COPY  (15*6) 

CALL  UNLOAHIIA) 

call  exit 

801  CONTINUE 

READ  INPUT  TAPE  5* 1 1 10*NC *NS* ( LS IK ) *K*1 *NC ) 

READ  INPUT  TAPE  5* 1 1 10* ( MB ( J ) * J*1 *NC ) 

NC1«NC4.1 

READ  INPUT  TAPE  5* 1 1 10* ( MB ( I ) *  I >NC1 *NS ) 

1000  F0RMAT(2I5*6(5XA5) ) 

1001  FORMAT! 12A6) 

1010  FORMATI 10X*6E10*5) 

1020  FORMAT! I5*A5*6(E5*0*5X) I 
1030  FORMAT! 7E1 0*5) 

1100  FORMAT !7E 10*5 ) 

1110  FORMATI14T5) 

1200  F0RMAT!215*4E12*4) 
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1900  format C 10F5.0) 

2213  FORMAT! lOH  ) 

DO  2  I«1»NS 
2  MC(  n«MB(  I  > 

CALL  SORTMCCMCtNSI 
DO  600  TJ«1»TJX 
00  6  1  «  !•  NC 
CP(I)  «  0.1 
CPLASTfl)  «  0.0 
6  CBLASTfM  ■  0.0 
RT«R«T 

IFCT-1600.0) 700* 700*701 

700  K»1 
GOTO  706 

701  IF! T-6500*0) 702*702*703 

702  K«2 

GO  TO  706 

703  IF! T-9500.0) 704*704*705 

704  K-3 

GO  TO  706 

705  K«4 

706  DO  22  l-l *NS 
M«MB! I  ) 

IF!ABSF!EKTH!K*H) ) -0*00001 > 707* 707  * 708 

707  WRITE  OUTPUT  TAPE  6*1201*T*M 
GO  TO  100 

700  EK0N—ATH!IC*M)*!  1*0-L0GF!T|  l♦EiCTH!K*M) 

TEMP«!  !  !  !  ETH!  K  *M  )  /4*0  )«T4>0TH!  K*M  )  /  3*0  )  •T^CTH !  K  *M  )  /2  *0  )  •T*>BTH!K  •M  I  I 
1*T4-EIC0N 
22  F2! I l«-TEMP 
r  ECU*  2 

00  28  ll«NCl*NS 
MI«MB!Iii 
DO  24  K1«1*NC 
LX«LS!IC1  ) 

DO  24  J1*1*NC 
MJbMB! J1 ) 

24  AIK1*J1I-ALF!LK*MJ) 

DO  26  1«1*NC 
LK-LS! 1 ) 

26  B! I laALFILKvMI ) 

CALL  MATIMV!A*NC*B*1 *OET ) 

IF!ABSF!0ETI-EPS2)30*30*32 
30  WRITE  OUTPUT  TAPE  6*3000*0ET 
GO  TO  lOO 
32  DO  28  J>1*MC 
28  ENU! J*I1 )«B! J) 

C  E0U*4 

DO  38  J«1*NC 
LJ>LS! J) 

DO  38  KK*1*NC 
MK«MB!KK I 

38  A! J*KKI«ALF!LJ*MK) 

DO  40  ICIC*1  *NC 
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LK«LSCICICI 

40  iicicici«rurficLKi 

CALL  NATINVCA»NC«B*l*OETt 
•FCABSFC0CT}-€PS2IA2*42«44 

42  MRITE  OUTPUT  TAPE  4«S100»0£T 
GO  TO  100 

44  DO  44  J«1«NC 
46  CUECJI«BfJI 
.C  EOU«9 

ENZ«0« 

DO  4B  KK«l«f«C 
L«LSCKKI 

4B  ENZ«EN24B£TAILI*CUEBCLI 
IFCSSW1.0.9}41«41«43 

43  WHITE  OUTPUT  TAP£6 •4000* C FZ I  I  I • I ■ 1 •NS  I 
DO  43  J«NC1«NS 

43  WHITE  OUTPUT  TAPE4 *4002 • C ENUI I • J I • I ■! •NCI 
WHITE  OUTPUT  TAPE4^4004^ ICUE f I  I • I ■! •NC I 
WHITE  OUTPUT  TAPE6^400S •ENZ 

41  CONTINUE 

hho«hhohz 

DO  300  JK«1  •  JICX 
00  36  laNCl^NS 
TEMP«0^ 

TMPB>0« 

OO  34  J«1*NC 

MaMnf  J| 

TENP«TEMP4ENUC J^I l•CHZ2CK•M|/HT♦FZC Jl I 
34  TMPB«TMPB^ENU( J^l I 

TNPB«ITMPB-1*>«L0GFI T/TBZ^HHOI 
M»MBf I  I 

TEMPaTEMP-HZZf  K^MI/HT-FZC  1  I^THPO 
34  EKfll«TEMP 

IFCSSW3-0^9 131 •31^33 

33  WHITE  OUTPUT  TAPE6^4004^ C EKI 1 1 • I aNCl •NS I 
31  CONTINUE 
300  NN  ■  1 

NCON  •  0 

C  BEGIN  LOOP  TO  CONVEHGE  INOEP*  SPECIES  FOH  GIVEN  T  AND  HHO/HHOZ 

-^Ol  DO  332  L  ■  !•  LIMIT 
fTEHfJKI  ■  L 

C  OOUBlE  LOOP  TO  CALC*  DEPEND*  CONCENTHAT IONS 

302  00  304  I  •  NCI*  NS 
TMP  •  ElC(  I  I 

OO  303  J  ■  1*  NC 

303  TMP  •  TMP  ♦  LOGFCCPC Jl I^ENUI J*l I 

304  cell  ■  EXPFl TMPI 

IF  CI4N  -  21  303*  303*  130 
303  OO  299  J  •  1*  NC 
24«j  KEVfji  •  3 
NOCON  •  0 

C  LOOP  TO  CALC*  CBIJI  AND  SET  CONVG*  LOGIC  AS  HEOO* 

DO  319  J  •  1*  NC 
TMP  .  0*0 


•3 


c 


306 


C 

c 


30T 


30S 

309 

310 

311 

312 

316 

3  16 

317 
3ia 
319 


320 


321 

606 

369 

350 

351 

352 

322 
C 


••23 


326 


607 

325 

611 

c 

326 


OBTAIN  CBfJI  USING  NOLE  OB  ELECTRON  BALANCE 

OO  306  I  «  NCI*  NS 

TMP  «  T^  ♦  ENUCJ*II6CCII 

SUMTIJ)  «  TMP 

CUEfJI  »  0*0  DENOTES  ELECTRON  BAL*  INO*  NOLES  NOT  FIXED  I 

CHFCX  NEG.  CUECJI 

IF  fCUECJM  307*  307*  308 

CBfJI  •  -  TMP 

CBCPIJI  •  CBCJI/CPCJI 

•CEVIJI  «  6 

GO  TO  319 

CBfJI  •  f fCUEf Jl /ENZ |6CPf Jl l/CCPf J I  ♦  TMP I 
CBCPfJI  ■  CBCJI/CPfJI 

IF  fABSFf CBCPfJI I  -  OCONVGl  312*  312*  311 
CfJl  •  CPfJI^OCONVG 
GO  TO  316 

IF  fABSFf CBCPf Jl  I  -  f 1 •O/OCONVGI I  313*  316*  316 
CfJI  ■  CPfJI/OCONVG 
AEvfji  •  1 


GO  TO  318 

IF  fL  -  II  319*  317*  319 
iCFVfj)  •  2 
NOCON  •  1 
CONTINUE 

LOOP  TO  CALC*  OOf^>  AND  CfJI  MMEN  DATA  IS  AVAIL* 

DO  322  J  •  1*  NC 

IF  fXEVfJ)  -  21  320*  321*  321 

QAfJi  •  0*0 

OOfJl  •  0*0 

GO  TO  322 


OAfJi  •  fCBfJI  -  CBLASTf Jl l/fCPI Jl  -  CPLASTfJil 

IF  fXEVfJ)  <-  61  369*  606*  369 

IF  fNCON  -  fNC^lll  351*  369*  351 

OQfJI  •  OAfJI/fGAfJI  -  1*01 

IF  fOOfJ)*  0*991  350*  351*  351 

IF  fOOfjll  351*  351*  352 

OOfJI  •  0*5 

CfJI  ■  OOfJIPCPfJI  6  f l*0-OOf Jl |6CBf Jl 
CONTINUE 


f <EVf Jl 


LOOP  TO  CALCULATE  CONVERGENCE  PARAMETERS 

00  323  J  •  1*  NC 

OELfJ)  ■  CfJI  -  CPfJI 

CONVGfJl  •  ANSFf DFLf Jl/Cf Jl I 

Sr.  IP  LOOP  TO  CMECIC  CONVERGENCE  IF  NOCON  •  1 

IF  fNOCONI  326*  326*  326 

NCON  •  0 

00  325  J  •  1*  NC 

IF  fCONVGfJI  -  EPSI  607*  325*  325 
NCON  •  NCON  ♦  1 

continue 

IF  fNCON  -  NCI  326*  611*  326 
MN  ■  5 

CHECIC  SSlf2  FOR  CONVERGENCE  TRAPS 
IF  f5SW2  -  0*51  610*  610*  327 
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127  WRITE  OUTPUT  TAPE  6*  5010*  L*  T*  RHO 

write  output  tape  6*  4011*  CSUMTfX).  I  »  1*  NO 
WRITE  OUTPUT  TAPE  6*  4012*  CCBtll*  I  •  1*  NCI 
WRITE  OUTPUT  TAPE  6*  4013*  fCBLAC^Tfll*  I  •  1*  NCI 
WRITE  OUTPUT  TAPE  6*  4014*  (CPCII*  I  •  1*  NCI 
WRITE  OUTPUT  TAPE  6t  4015*  fCPLASTIII*  I  •  1.  NCI 
WRITE  OUTPUT  TAPE  6*  4016*  (OAfll*  I  «  1*  NCI 
WRITE  OUTPUT  TAPE  6*  4017.  fOOtll*  I  •  !•  NCI 
WRITE  OUTPUT  TAPE  6*  5011*  (ICEVdl*  I  ■  NCI 
WRITE  OUTPUT  TAPE  6«  4018*  fCBCPfll*  I  •  1*  NCI 
WRITE  OUTPUT  TAPE  6*  4019*  (OELfll*  I  ■  1*  NCI 
WRITE  OilTPUT  TAPE  6«  4020*  lOIl*  I  ■  1*  NCI 
WRITE  OUTPUT  TAPE  6*  4021*  CCONVGf 11 •  I  «  1*  NCI 


C  SKIP  TRIAL  VALUES  IE  LOOP  DOES  NOT  CONVERGE 

410  IE  (L  -  LIMITI  328t  332.  332 

C  SET-UP  trial  values  EOR  next  ITERATION  OR  INITIAL  VALUES  NEXT  RHO 

328  00  331  J  »  !•  NC 

IE  fKEYlJI  -  21  330.  329.  330 

329  CPLASTCJI  ■  CP(JI 
CPiji  ■  rpfji 
CBLASTIJI  «  CBfJI 
GO  TO  ‘•3  1 

330  CPLASTCJI  «  CPCJI 
CPCJI  ■  CCJI 
CRLASTCJI  «  CACJl 

331  continue 

C  IE  CONVERGED. CNN  •  91  RETURN  TO  302  AND  CALC.  DEPEND.  CONC.  ONLY 

IE  (WN  -  21  332.  332.  302 
C  END  OE  CONVERGENCE  LOOP 

332  CONTINUE 

C  IE  limit  exceeded.  PRINT  CURRENT  RESULTS 


WRITE 

WRITE 

OUTPUT 

OUTPUT 

TAPE 

TAPE 

6  . 
6. 

3300. 

4011. 

T.  RHO 

C  SUMT  f I  1  . 

I 

m 

1 . 

NCI 

WRITE 

OUTPUT 

TAPE 

6  . 

4012. 

CCBCII. 

I 

m 

1 . 

NCI 

WRITE 

OUTPUT 

tape 

6  . 

4013. 

CCBLASTf I  1  . 

I 

m 

1 . 

NCI 

WRITE 

OUTPUT 

TAPE 

6. 

4014. 

C  CP  f I  1 . 

I 

m 

1 . 

NCI 

WRITE 

OUTPUT 

TAPE 

6  . 

4019. 

(CPLASTC I  1 . 

I 

m 

1 . 

NCI 

WRITE 

OUTPUT 

TAPE 

6  . 

4016. 

fOAf I  1 • 

I 

m 

1 . 

NCI 

WRITE 

OUTPUT 

TAPE 

6  . 

4017. 

fOOC I  1 . 

I 

m 

1. 

NCI 

WRITE 

OUTPUT 

TAPE 

6  . 

5011. 

CKEYI I  1 . 

I 

m 

1 . 

NCI 

WRITE 

OUTPUT 

TAPE 

6  . 

4018. 

CCBCPf I  1 . 

I 

m 

1 . 

NCI 

WRITE 

OUTPUT 

TAPE 

6  . 

4019. 

COELC I  1 . 

I 

m 

1 . 

NCI 

WRITE 

OUTPUT 

TAPE 

6  . 

4020. 

(Cf  1  1. 

I 

m 

1. 

NCI 

WRITE 

OUTPUT 

TAPE 

6  . 

4021. 

CCONVGf  I  1  . 

I 

m 

1 . 

NCI 

C  ZERO  OUT  OUTPUT  ARRAY 

GO  TO  “^AO 
130  Z  •  0.0 

no  132  J«1.NS 
132  2  ■  2  ♦  CCJI 

P  ■  2«RM03C T/TB2 I 
ENV*RH03Z3ENL 
DO  134  J-l.NS 
PCCCJI  ■  CCJ|3RH0*ENL 
134  XCJ|«CfJI/Z 
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STORE  IN  OUTPUT  ARR4.V 
J3«JK4>1 
OAf J3»l |«RHO 
OPf J3«l |«RHO 
OAf J3»2|«P 
OPC J3*?|»P 
OAC J3»3l«Z 
OPfJ3»3l»Z 
OA( J3*4)«ENV 
OPfJ3«4l«ENV 
DO  136  J«1*NS 
J1«MC( Jl 
J2-MRIJ1  14^4 
OPCJ3«J2)  *  PCCIJll 
136  OAfJ3tJ2) «X( Jl ) 

500  RHO«RHO«DRH 
560  M*JICX4.1 

L«6 

WRITE  OUTPUT  TAPE  6*2207 

501  WRITE  OUTPUT  TAPE  L*2213 
WRITE  OUTPUT  TAPE  L*2200*T 

WRITE  OUTPUT  TAPE  L*2202 • C 0A( I • 1 ) • I «2 tM ) 
WRITE  OUTPUT  TAPE  L • 2204 » < OA f I *2 ) • I *2 ) 
WRITE  OUTPUT  TAPE  L* 2206* I OA ( I • 3 ) • I -2 tM ) 
WRITE  OUTPUT  TAPE  L • 2206 * ( OA ( I *4 ) • I *2 tM ) 
WRITE  OUTPUT  TAPE  L*2212 
WRITE  OUTPUT  TAPE  L*2209 
IF(K-2I503*507*507 

503  DO  509  J»1*NSP 
J1«J4.4 

IF(OA(2* Jl ) 1 509 *505 *504 

504  WRITE  OUTPUT  TAPE  L*2210* (OAf I • Jl ) *  I«1  *M  ) 

505  CONTINUE 
M«JKX^1 
L«6 

K«1 

WRITE  OUTPUT  TAPE  6*22<>7 


0 

6  1 


E  O 


UTpUT  T 


aPe  *-•2213 


WRITE  OUTPUT  TAPE  L*2200*T 
WRITE  OUTPUT  TAPE  L * 2202 • f OP ( I • 1 ) *  I «2 vM ) 
WRITE  OUTPUT  TAPE  L * 2204* f OP f I *2 ) • I *2 I 
WRITE  OUTPUT  TAPE  L • 2206* f OP f I • 3  I  *  I *2 tM ) 
WRITE  OUTPUT  TAPE  L*2208 • I OP f I *4 ) • I *2 I 
WRITE  OUTPUT  TAPE  L*  2214 
WRITE  OUTPUT  TAPE  L*2209 
IFIK-2»603*907*507 
603  DO  605  J«1*NSP 
J1-J4.4 

IFCOPf 2* Jl I ) 605 *605 *604 
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604  WRITE  OUTPUT  TAPE  L • 22 10* COP C I • J1 I • I •! I 

605  continue 

WRITE  OUTPUT  TAPE  6*  4023*  CITERCll*  I  •  1*  JKX) 

506  L«15 
K*5 

WRITE  OUTPUT  TAPE  15«2211tJICX 
60  ro  501 

507  DO  508  J«1»NSP 
J1»J4>4 

508  WRITE  OUTPUT  TAPE  L*2210* COAI I •Jl I • .N) 

DO  799  J«1»NSP 

DO  798  I  ■  2*  M 
OPC I • Jl  ■  0«0 
798  OAfIvJ)  >  0*0 
799  CONTINUE 
600  T»T^OT 

GO  TO  100 

1201  FORMAT (43HOTHERMOOYNAMIC  FUNCTION  NOT  AVAILABLE  T  -1PE13«4«13H 
1  SPECIES  *131 
2200  FORMATf6HO  T  «F6«0) 

2202  FORHATCOHO  RHO  •(lPllEll«3n 

2204  FORMAT(8HO  P  •llPllEll«3n 

2206  FORMATCOH  2  • C IPl lEl 1* 3 ) I 

2207  FORMAT! INI) 

2208  FORMATC8M  NV  • f IPl lEl !• 3  I  I 

2209  FORMAT! IHO) 

2210  FORMAT! 1A6«2X«! 1P11E11«3 n 

2211  F0RMAT!lHltl4) 

2212  FORMAT C45H0SPEC IE  CONCENTRATIONS  IN  MOLE  FRACT lONS-P ! S M 
2214  FORMAT C48H0SPEC IE  CONCENTRATIONS  IN  PART ICLES/CM4«3-PCC ! S ) I 
3000  FORMAT! 1HO/28HOEOU.2  FAILURE  DETERMINANT* tlPElI .A) 

3100  FORMAT  ! 1HO/28HOEOU«4  FAILURE  DETERMINANT- • IPEl 1 *4  ) 

3300  FORMAT!21HONO  CONVERGENCE  T  -IPElS.SvSH  RHO  -1PE13.6) 

4000  FORMAT!7hOFZ!N5I  /! IH  *iP8E 14«5 1) 

4002  FORMATC4H  ENU  / !  IH  1P0E14«5M 
4004  format  !8H0EIC!  NS- )  /!1H  1P8E14.5II 
4006  FORMAT !8H0CUE! NCI  /!1M  1P8E14.5II 
4008  FORMAT!5HOENZ«1PF12«5I 

4010  FORMAT! 1X«1P8F14«5I 

4011  FORMATC7M  SUMT  •C1P8E14.5II 

4012  FORMAT! 7H  CB  •!1P8E14«5II 

4013  FORMAT! 7H  CBLAST  •C1PBE14.5I1 

4014  FORMATC7H  CP  •!1PBE14«5II 

4015  FORMATC7H  CPLAST  •!1P8E14«5II 

4016  FORMAT! 7H  OA  •!1P0E14«5II 

4017  FORMATC7H  00  •!1P8E14«5II 

4018  FORMAT! 7H  CBCP  •C1P8E14.5II 

4019  FORMAT!  7H  DEL  •C1P8E14.5M 

4020  FORMATC7M  C  •!1P8E14«5II 

4021  FORMAT! 7M  CONVG  •CIPBEIA.SM 
4023  FORMATC6HO  ITER*2X»  C9111II 

5010  FORMAT! 1H0*9X»3HL  •  I4»7X«3HT  •  0PF8*1«4H  R-  1PE10«3I 

5011  FORMATC7N  XEV  •C8I14II 
END 
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matrix  INVCKSION  WITH  ACCOMPANYING  SOLUTION  OF  LINEAR  EQUATIONS 
SUBROUTINE  MAT INVC A*N»B«M»OETERM) 

DIMENSION  IPfVOT f 10 1 •A C 10* 10 1 •Bf 10 • 1 ) • INDEX C 10*2  I •PIVOTC 10| 
equivalence  flROWsJROWI*  C ICOLUM* JCOLUM ) •  (AMAX*  T*  SWAPi 

initialization 

10  0ETERM«1«0 
IS  DO  20  J«1*N 
20  1PIV0TCJ)«0 
30  DO  SSO  I«1»N 

c 

C  SEARCH  FOR  PIVOT  ELEMENT 

C 

40  AMAX«0«0 
4S  DO  lOS  J«1«N 

SO  IF  (  IPiVOTf Jl-1 I  60*  10S«  60 
60  00  I on  X«1*N 

70  IF  ( 1P|V0T(KI-1 I  80*  100»  740 

80  IF  f ABSFf AMAX|-ABSF( A( J«KI 1 )  8S»  100.  100 

8S  IROW-J 

90  ICOLUMaK 

9S  AMAX«AfJ«K) 

100  CONTINUE 
lOS  CONTINUE 

110  IPIVOTI ICOLUMIalPlVOTC ICOLUMi^l 
r 

C  INTERCHANGE  ROWS  TO  PUT  PIVOT  ELEMENT  ON  DIAGONAL 

130  IF  ( IROW-ICOLUM)  140*  260»  140 
140  DETERMa-DETERM 
ISO  00  200  L«1*N 
160  SWAPaAC IROWtLl 
170  AC IROW*L)-A( ICOLUMtLl 
200  A( ICOLUMvLIaSWAP 
20S  IF(MI  260*  260t  210 
210  00  2S0  L«lt  M 
220  SWAPaBf IROWvLI 
230  B( IROWtL I -Bf ICOLUMvL) 

2S0  B( TCOLUM«L)aSWAP 
260  INOFXC I •! lalRow 
270  INDEX! I *2 l>ICOLUM 
310  PI  VOT!  DaA!  ICOLUMt  ICOLUM) 

320  OETERMaOETERM^PlVOT!  n 

01 VI OF  PIVOT  ROW  BY  PIVOT  ELEMENT 

330  A! ICOLUM* ICOLUM) al«o 
340  DO  3S0  L«1«N 

3S0  AC ICOLUM*L)aA( ICOLUMtL )/PlVOTf I ) 

3SS  IF(M)  380*  380*  360 
360  DO  370  L«1«M 
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970  BC ICOLUM«LI«BC lCOLUM»L I /PI VOT C I  I 
C 

C  REDUCE  NON-PIVOT  ROWS 

C 

380  DO  950  L1«1*N 

390  IFCLl-ICOLUMI  400*  950*  400 

400  T«A(L1«IC0L0M) 

420  ACL1*ICOLUMI«0«0 
430  DO  450  L«1»N 

450  AlLl«LlaAlLi«Ll-Af IC0LUM*L)*T 
455  IF(M}  990*  590*  460 
460  DO  500  L«1«M 

560  Bf Ll*L)«VrLl*L)-B( ICOLUNtL)*r 
950  CONTINUE 
C 

C  INTERCHANGE  COLUMNS 
C 

600  DO  710  laltN 
610  L«mi-T'  ■ 

620  IF  ( IN0EX(L«ll-lNDEXfL*2M  630*  710*  630 
630  JROWalNDEXIL*!} 

“  640  JCOUuHiTNOEXTLii I  ~ 

690  DO  705  K«1«N 

660  SWAPaACKtJROWI 

6T0  A(IC» jllOW)«ACK«JCOLUM| 

700  A(K«JCOLUM)aSWAP 
709  CONTINUE 

~T10  WNTTNUf - 

740  RETURN 
END 

SUBROUTINE  SORTNECMO^S) 

DIMENSION  MCCll*M0f50l 
DO  2  laltNS 

2  NDM)-! - - - 

Kl-NS-1 

DO  2«  JIalvXl 

t2«Jl4l 

DO  20  ll«l2tNS 

IFCMCC J1 l-MCC I1II20»20«10 

10  ir-Mcurr  ■ 

MCCJlIaMCdll 
II  l^£T 
IT2-MDCJ1  I 
MD( JllaMDCIll 
MOC  IH-IT2 
20  CONTINUE 

DO  30  lalvNS 
30  MCdtaMDf  II 
RETURN 
END 
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